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STAGNATION- POINT SOLUTIONS FOR INVISCID 
RADIATING SHOCK LAYERS 1 


By Walter B. Olstad 
Langley Research Center 

SUMMARY 

Four approximate analytic solutions are presented which are valid in the following 
situations: (1) the effects of radiation can be treated as small perturbations; (2) the opti- 
cal thickness of the shock layer is small; (3) the optical thickness of the shock layer is 
large; and (4) radiation is the principal mode of energy transport. The first two solutions 
include the effects of variable thermodynamic properties, variable and nongray optical 
properties, and body surface reflectivity. The third and fourth solutions are restricted to 
gray gases. 

The various solutions were used to indicate the effect of the radiation cooling param- 
eter, the Bouguer number, the enthalpy dependence of the absorption coefficient, and sur- 
face reflectivity on the shock-layer enthalpy distribution, the stagnation-point radiant heat- 
transfer rate, and the shock standoff distance. A discussion of the effect of radiation 
cooling on the stagnation-point convective heating rate is also included. The results indi- 
cate that at every altitude and velocity, there is a finite value of body nose radius for which 
the rate of radiant heat transfer is a maximum. Also a significant reduction in the com- 
puted value of the radiant heating results when the nongray character of the absorption 
coefficient is taken into account. 


INTRODUCTION 

As the exploration of space progresses from the near-earth environment to the moon 
and the planets of the solar system, study of the atmospheric entry of objects in excess of 
escape velocity (about 11 km/sec) becomes necessary. In addition to studies concerning 
man-made objects, there is considerable interest in the entry of meteoroids into the 


i-This report is a slightly revised version of a thesis presented to the Division of 
Engineering and Applied Physics, Harvard University, in partial fulfillment of the require- 
ments for the doctor of philosophy degree in Applied Mathematics, April 1966. 



earth's atmosphere at velocities from 20 to 70 km/sec. At these large speeds, radiant- 
energy transfer is an important factor governing the behavior of the hot- shock- layer gas 
enveloping the object. 

Consequently, a number of investigators have addressed themselves to the problem 
of the radiating shock layer. The first analyses assumed that the flow processes were 
unaffected by (or uncoupled from) the transfer of energy by radiation. (See, for example, 
refs. 1 and 2.) Thermodynamic and flow properties were calculated by neglecting radia- 
tion. The radiant energy flux was then calculated from measured or theoretically deter- 
mined optical properties for these conditions. Although this approach provided accepta- 
ble engineering estimates at speeds less than escape velocity, it was not sufficient to 
describe the effects of radiation at higher speeds. The next step was to take into account 
the loss of energy from the shock layer due to radiation. This cooling of the shock layer 
tends to reduce the emergent radiant energy flux. This reduction is often termed "radia- 
tion decay." Radiation decay was studied by a number of authors. (See, for example, 
refs. 3 to 8.) All the cited works, with the exception of reference 8, used the transparent 
approximation, 2 which neglects absorption within the shock layer, and assumed that the 
radiation cooling of the shock layer gas was small and did not influence the mass trans- 
port. The process of absorption by a gray 2 gas was studied in reference 8. However, 
the flow model used in that investigation only roughly approximates the flow in the stagna- 
tion region of a shock layer. Consequently, the analysis was unable to describe details of 
conditions in the shock layer or to provide reliable quantitative results. 

Perhaps the most ambitious analysis to appear to date is the work of Howe and 
Viegas (ref. 9). (This statement was true at the time of first writing. Since then, how- 
ever, a number of papers which treat the same complexities, some in more detail, and 
which include nongray effects have been published.) They obtained numerical solutions 
to the integro- differential system of equations governing the flow in the stagnation region 
including the effects of radiation decay, absorption by a gray gas, viscosity, and surface 
mass injection. 

All the works discussed are restricted to velocities less than about 20 km/sec, 
although the work of Howe and Viegas was so restricted simply because they did not 
choose to make calculations for higher velocities. Fay, Moffatt, and Probstein (ref. 10) 
undertook an analysis of meteoroid entry, in the speed range of 20 to 70 km/sec. Since 
they were interested only in obtaining upper bound estimates of radiant heating, they 


2it was so-called because the shock- layer gas is considered to be transparent to its 
own radiation. 

3 A gray gas is one for which the optical properties are independent of the photon 
energy or wavelength. 
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ignored radiation decay and absorption (except that they did not allow the radiant energy 
flux to exceed the black-body limit), both of which can be very important at these speeds. 

Although the existing studies (which include many works in addition to those cited) 
have contributed much to the qualitative and quantitative understanding of the physical 
processes taking place in radiating shock layers, a great amount of work remains. For 
example, parametric studies of absorption in a realistic shock-layer flow are lacking, the 
effects of surface reflectivity have been generally ignored, and there have been no reported 
attempts (at least in the knowledge of this investigator) to study shock-layer gases with 
nongray optical properties. At the time this paper was first prepared Lick (ref. 11) and 
Greif (ref. 12) considered nongray optical properties in their studies of combined radiation 
and conduction. Their results indicate that nongray effects can be significant. 

The investigation reported herein was undertaken to provide a parametric study of 
the influence of radiation on blunt objects large and small, traveling at speeds up to 
70 km/sec. The approach was to seek simple approximate solutions where available in 
the hope that they would lead to a better understanding of the physical processes involved. 
The parameters studied included the radiation cooling parameter e (which characterizes 
the relative importance of radiation as an energy transport mechanism compared with 
convection), the Bouguer number (which indicates the importance of absorption in radiant 
transport), the surface reflectivity (which is indicative of the ability of the surface of the 
object to accept the incident radiant energy flux), and the spectral variation of the absorp- 
tion coefficient. (There is no single quantity or even group of quantities which character- 
izes the important effect of spectral variation on the flow.) Definitions of these parame- 
ters and their role in influencing the flow is discussed in greater detail in subsequent 
sections. 

In order to facilitate this investigation without sacrificing physical significance, the 
analysis was limited to the stagnation region and the following conditions were assumed 
to apply: (1) the shock-layer gas is in local thermodynamic and chemical equilibrium, 

(2) the body geometry is axisymmetric, (3) there is no mass addition to the flow from the 
body surface, (4) the thicknesses of the shock and the viscous boundary layer are small in 
comparison with the shock_ standoff distance, and (5) absorption in the free stream ahead 
of the object is negligible. 

In this investigation, solutions were obtained for four limiting cases of the radiation 
cooling parameter and the Bouguer number. The first of these is a small-perturbation 
expansion in the radiation cooling parameter e, which is valid when the influence of radia- 
tion is small. The second solution holds when the shock layer is optically thin. This 
solution is presented as a small-perturbation expansion in the Bouguer number. A solu- 
tion valid when the shock layer is optically thick (Bouguer number » 1) and the final 
solution, which is restricted to the case when radiation is the principal mode of energy 
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transport within the shock layer, are also presented. The first and second solutions have 
been formulated to include the effects of nongray radiation. The third and fourth solu- 
tions are restricted to the gray case. In each of the four limiting cases, it is possible to 
approximate the governing integro-differential system of equations by a purely differen- 
tial system which leads to a singular perturbation problem. 

The results obtained by means of the various approximations are combined to give 
the radiant-heat-transfer rate and an estimate of the effect of radiation on the convective 
heating rate at the stagnation points of blunt objects traversing a gray model earth atmo- 
sphere. The effects of the nongray character of air on these results are discussed. 

In appendix A, a boundary -layer analysis is performed on the integro-differential 
system to determine the form of the equations in the inviscid region and the viscous 
boundary layer. Mathematical developments for optically thin shock layers and radiation- 
depleted shock layers are presented in appendixes B and C. The symbols used in the 
analyses are given in appendix D. 

Since this work was completed in early 1966 a number of studies which include the 
effects of radiation decay and nongray absorption have appeared in the literature. The 
importance of atomic line radiation was recognized and has been included. The influence 
of surface ablation and absorption in the free stream have also been studied. Excellent 
discussions of these more recent developments can be found in two survey papers 
(refs. 13 and 14). The author is indebted to Professor H. W. Emmons for his guidance 
during the course of this study and to Mrs. J. T. Kemper who did most of the computer 
programming. 

STAGNATION MODEL FOR A RADIATING SHOCK LAYER 

Fundamental Equations of Radiation Gas Dynamics 

Prior to setting up a particular flow model for the problem at hand, it is desirable 
to examine briefly the fundamental equations of radiation gas dynamics. An excellent 
discussion of these equations has been presented by Goulard (ref. 15), and the reader is 
referred to this work for a more detailed exposition. 

The studies of this paper are limited to the steady flow of gases in local thermody- 
namic and chemical equilibrium. In addition, the effects of radiation pressure and radia- 
tion energy density are ignored. These effects are important only when the radiant energy 
flux is extremely large as it is deep in the interior of a stellar atmosphere. Finally, the 
presence of external forces, such as gravity and electromagnetic forces, are neglected. 
With these restrictions in mind, the conservation equations for a radiating gas can be 
written 
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H,i = 0 


(Continuity) 


( 1 ) 



(Momentum) 


( 2 ) 


vR (Energy) 

where the quantity h^ is the total specific enthalpy of the gas 

h t = h+ i u i u i 


^j u u = p,i + T ij,i 


pU i\i = -( U i T ii),j ' q i,i " q i,i 


(3) 

(4) 


The static specific enthalpy h includes the chemical energy of the gas in terms of the 
heats of formation of the various gaseous species. 

The double subscript notation is employed wherein a subscript (i,j,k) denotes the 
(ith, jth, kth) component of a vector, a subscript preceded by a comma denotes differen- 
tiation of the subscripted quantity with respect to the (ith, jth, kth) direction, and any 
term which displays a repeated subscript represents the sum of terms over all possible 
values of the repeated subscript. 

An expression relating the thermodynamic variables is needed to complete the set 
of equations. A convenient form is 


h = h(p ,p) 

The molecular transfer processes are represented by the classical expressions 

T ij = M ( u u +u j,i) + ( ,J '-t ( 1 ) 6 ij u k,k 


(5) 

( 6 ) 


= - k eff T ,i 


(7) 


The quantity k e jf is an effective coefficient of heat conduction which includes the effects 
of energy transport by molecular collisions and by the diffusion of reacting species. 

These two processes can be lumped together like this only when the conditions of local 
thermodynamic and chemical equilibrium hold,_ (See ref. 16.) 

D 

The radiant-energy-flux vector qj is defined in terms of the radiation inten- 
sity J x 









q 


R 

Xi 



dw 

J 


( 8 ) 
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and is the rate of flow of radiant energy per unit area across an element of area whose 
normal points in the ith direction. The quantity l j is the direction cosine between the 
direction of a single beam of intensity and the ith direction. Then J x can be deter- 
mined from the conservation equation of radiation transfer 


-W=-P^X-f x ) (9 

where /3 x is the mass extinction coefficient. It is composed of the mass absorption 
coefficient k x and the mass scattering coefficient a x 

do: 

The ratio of mass emission coefficient j x to the mass extinction coefficient j3 x is 
often called the source function S x = 

For n onsc atte ring me dia-ina state of local thermodynamic equilibrium (ct x = o), the 
source function reduces to the Planck function 

B x = 2hc2^hc/XkT _ j) 1 (H) 


provided that the mass absorption coefficient k ^ includes the effects of induced emis- 
sion. Here h and k are the Planck and Boltzmann constants, respectively, and c is 
the speed of light. Throughout the remainder of this paper, it is assumed that the gas in 
the shock layer is nonscattering. This assumption is reasonable as the number of large 
solid particles which might scatter radiation is expected to be negligible in the shock 
layer. A few such particles might exist in the cooler regions of the boundary layer 
adjacent to the body surface as a result of "spalling" of this surface. However, their 
presence might be accounted for, if necessary, by changing the effective reflectivity of 
the body surface. 

At the extremely high shock-layer temperatures for which the gas is multiple ion- 
ized and free electrons are plentiful, Thomson scattering can become important. For 
example, Kivel and Mayer (ref. 17) show that scattering cannot be neglected when the tem- 
perature reaches 350 000° K at densities less than about 0.01 of the sea-level value. 

For the nonscattering case the intensity of radiation at a point M in the direc- 
tion s follows from a formal integration of equation (9) 


fs) (s) 

h = J p=Q B x( p ) ex P( -T x,Mp) dr X,MP + J X exp (“ r X,MQ) 
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where 


r X,MP " 



ds 


(13) 


is the "optical thickness" or "optical path length" along the beam between the points M 
and P and P is a "running" point on the beam between point M and the boundary point Q. 
(See sketch (a).) Although the terms optical thickness and optical path length, long 
established in astrophysical literature, seem to imply a dimension of length, the quan- 
tity jyjp is dimensionless and is indicative of the number of photon mean free path 
lengths in the physical distance between M and P. 



Sketch (a) 

The quantity j| S \q) represents the contribution to the intensity at point M from 
the boundary and, in general, includes emission from the surface, reflection from the 
boundary of radiation originating from within the region, and transmission through the 
boundary of radiation originating from without the region. 

The integral term represents the summation of the contributions from all points P 
along the beam reduced by the attenuating factor exp^-T^ which accounts for absorp- 
tion by the intervening matter. 

The divergence of the radiation flux vector can be found with the aid of solution (12) 
with the result 


<1 - - JTMm J v K b x)p ^pM dv(p > * 


So K)mI a j ‘ s><q) 


( s ) //->\ exp (~ r X ,Mq) 


MP 

cos 9 ds(Q) dX 


(14) 


MQ 
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The integrations over the volume V of the gas and the area A of the bounding surface 
include only those portions of the volume and surface which are visible to an observer 
stationed at point M. 


Stagnation Flow Model 

The study of three-dimensional flow of compressible gas in the vicinity of the for- 
ward face of a blunt body cannot be reduced by transformation to the study of an equiva- 
lent one -dimensional system as can be done in the incompressible case. However, avail- 
able numerical solutions (see, for example, refs. 18 to 20) indicate that for all practical 
purposes, a reduction from a three-dimensional to a nearly equivalent one-dimensional 
problem can be carried out in the stagnation region. The reason that this simplification 
can be applied is thatthe coupling between the momentum and energy equations is very 
weak. Also the various thermodynamic properties are nearly independent of the lateral 
or radial coordinate. Although the same arguments apply in the stagnation region of a 
radiating shock layer, it is not possible to postulate the existence (even approximately) of 
a one -dimensional solution solely on this basis. Some additional assumption is required 
regarding the effect of the far field on the radiant heat flux and its divergence. This 
effect, of course, cannot be obtained "a priori" as it depends on the solution to the entire 
flow field. Fortunately, the shock layer is thin and only a small part of the radiant energy 
emitted by gas removed from the stagnation region actually passes through the stagnation 
region. If absorption is small, only a small part of this energy is absorbed in the stagna- 
tion region. If, on the other hand, absorption is large, the beam is greatly attenuated when 
it reaches the stagnation region and thus only a small part of the energy which started the 
journey remains to be absorbed in the stagnation region. The divergence of the radiant 
flux is influenced only by the amount of energy absorbed and emitted. Consequently, the 
far-field effect on the divergence of the radiant flux is a result of that small part of radi- 
ant energy originating in the far field and absorbed in the stagnation region. In the trans- 
parent and optically thick limits, this effect of the far field vanishes. 

It would appear from this discussion that a stagnation model for a radiating gas can 
be postulated as long as the assumptions concerning the far field are not grossly unrealis- 
tic. In what follows, a particular stagnation model is formulated and an estimate of the 
inaccuracy resulting from the assumption concerning the far field is obtained. 

A schematic of the flow in the stagnation region of a blunt body is shown in sketch (b) 
At very high speeds, the ratio of the shock standoff distance A to the shock radius R s 
is very much smaller than one (a typical value is 0.05). Under these conditions, the geom- 
etry of the stagnation region closely resembles a plane parallel gas slab. In addition, the 
enthalpy in the shock layer varies slowly with respect to r/A so that the stagnation 
region may be approximately represented by a gas slab in which the thermodynamics as 
well as the geometry is one dimensional. 
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Sketch (b) 

As a result of the above considerations, the model described below has been chosen 
to represent the flow of a radiating gas in the stagnation region of a blunt object. The 
model consists of an axially symmetric flow impinging upon an infinite flat plate normal 
to the stream direction. At a plane which is parallel to the plate and a distance A in 
front of it, the gas is suddenly raised to a total specific enthalpy of W^. The plate is 
held at a constant temperature T w . Sketch (c) illustrates the geometry of the flow model. 


u 



^,s = | W ' 

Shock 


r 



k 

/ 

/ 

/ 

/ 


/ 


/ 

T w 

/ 

w 

/ 


/ 


/ 


Body 


Axis of 
symmetry 


Sketch (c) 
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The general equations of motion (eqs. (1) and (3)) when specialized to the axisym- 
metric geometry become 


fr (pur) + fe (pwr) = 0 


9u . 8u _ 8p . 8r rr . 9r rz . T rr " T 66 

— — + pW — — — H 1 H 

8r 0 r a z r 


pu ^ + pw n 
^ 8r H 8z 


pU tF + pW 9z 


M + oW M=_M + t zz + ° T rz Irz 
p 8z 8z 8r r 


(15) 

( 16 ) 
(17) 


pu 


£^t + pW £iH := _^r _^r + 

9r p 8z 8r r 8z 8r\ 


8h t _ 8 q r <lr d( *z 8 


UT rr + W T rz 




Ur 4 . wt 

'rr r w 'rz 


HkC 


UT rz + WT zz 


) 

(18) 

where q r and q z are the r- and z-components, respectively, of the heat flux vector 
which includes conduction, diffusion of reacting species, and radiation. The stress com- 
ponents are given by the expressions: 


The equation of state is 


'rr -»«* &♦(!•' -f (*)(!? + If *S) 
= If 4' -!")(§?+ 
r„ - 2* S4 -f ■.)(&»♦!?♦») 

Tzr = + |u’ 

zr ^\8r 8z, 

h s h(p,p) 


(19) 

(20) 
( 21 ) 
( 22 ) 

, (23) 


In order to specify the problem completely, a consistent set of boundary conditions 
must be provided. The kinematical conditions on the velocity are 


w(r,0) = 0 


(24) 


p(r,A) w(r,A) = -p^ (25) 

The first of these conditions restricts the analysis to one for which there is no 
injection from the surface of the object. When gas injection is important, it is necessary 
to replace the zero on the right-hand side of equation (24) with w w , the normal velocity 
of the gas at the wall. The second condition was obtained from continuity across the 
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normal shock at the stagnation point. A third kinematical condition is introduced here in 
order to relate the standoff distance A to the variation of the tangential velocity along 
the surface z = A. This variation in velocity is taken to be equal to that behind the near- 
normal part of a spherical shock; that is, 

u(r,A) = Woo cos /3 ~ 

I 

'! where /3 is the local inclination of the shock from the free- stream direction and 

j the radius of the spherical shock. 

j The dynamical "no- slip” condition at the surface is 

i 

! u(r,0) = 0 

The conditions on the enthalpy and pressure are 

h t (r,A)=|wJ, (28) 

h t (r,0) = h w (29) 


(26) 
R s is 


(27) 


p(r,A) 




(30) 


where x - P fp( 0,A) is the density ratio across the normal shock. Condition (28) comes 
from the conservation of energy across a strong normal shock and does not take into 
account absorption in the free stream of radiant energy emitted by the shock layer. Con- 
dition (29) restricts the analysis to those conditions at which a temperature "jump" or 
discontinuity is not present at the body surface. Such a "jump" can occur only when the 
molecular mean free path in the gas is not negligible in comparison with the characteris- 
tic length of the domain (in this case, the thickness of the thermal boundary layer). Con- 
dition (30) is the pressure distribution behind the near-normal part of a spherical shock of 
radius R s . 

j In addition to the boundary conditions listed, boundary conditions on the radiant 

j energy flux must be specified. These conditions are: 

I (1) The boundary at z = A (which corresponds to a bow shock) is transparent 

I 

(2) There is no radiant energy transfer from the free stream to the shock layer 

(3) The boundary at z = 0 (which corresponds to the body surface) is cold and 
reflects diffusely and independently of wavelength a fraction r w of the incident radiation. 

The statement (contained in condition (3)) that the body surface is cold means that 
emission from the body surface has a negligible influence on the gas in the shock layer. 
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When the hot (temperatures in excess of 10 000° K) shock layer is optically thin, emission 
from the relatively cool (temperatures less than 4000° K) body surface may be comparable 
with emission from the shock-layer gas. However, because the shock layer is optically 
thin, very little of the radiant energy emitted at the body surface will be absorbed by the 
shock-layer gas. On the other hand, when absorption in the shock layer is important, the 
shock-layer gas emission will approach the black-body value corresponding to the high 
shock-layer temperature. Since black-body radiation is proportional to the fourth power 
of temperature, the gas emission from an optically thick layer will greatly exceed the 
emission from the body surface. Thus, whenever the body surface temperatures are 
small compared with the shock-layer gas temperatures, the influence of emission from 
the body surface on the shock-layer gas is unimportant. 

It can be seen from the definition of the total enthalpy 

fy. = h + i(u 2 + w^) 


and the boundary conditions (25) and (28) that the magnitude of the kinetic energy in the 
shock layer is of order y 2 compared with the static specific enthalpy. For a strong 
shock, which is the only case of interest here, 0.05 is a typical value for y, the density 
ratio across the shock. As a consequence of the small value of y, the kinetic energy 
terms are neglected in the subsequent analysis. The viscous dissipation terms (the last 
three terms on the right-hand side of eq. (18)) are also neglected because only kinetic 
energy is dissipated through the action of the viscous forces. 


It is desired that the solutions to the one -dimensional model represent, as closely 
as possible, the phenomena in the stagnation region of a blunt body. For simplicity, the 
blunt-body geometry, flow field, and thermodynamic properties are considered to be 
axially symmetric about the stagnation streamline. Expanding the solutions for the real 
blunt-body problem in terms of the radial coordinate r and arguing on physical grounds 
that w(r,z), p(r,z), p(r,z), and h(r,z) are even functions of r whereas u(r,z) is 
odd, gives 


w = w^(z) + 0(r 2 ) 
u = ru^(z) + 0(r 2 ) 

P = P^(z) + 0(r2) \ 
p = z) + 0(r 2 ) 
h = h^(z) + 0(r2) 


(31) 
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In addition, the heat flux components will have the form 


q z = q i°^( z ) +°(r 2 ) 
q r = rq^(z) + 0(r3) 


( 32 ) 


Neglecting terms of order r^ and higher restricts the solutions to the vicinity of 
the stagnation point. Since stagnation region solutions are desired, it is assumed that the 
solutions in the plane parallel model have the functional forms of equations (31) and (32) 
truncated after the linear term in r. For these assumed forms, the continuity equa- 
tion (15) requires 


pu = rg'(z) 
pw = -2g(z) 


(33) 


That portion of the heat flux due to conduction and diffusion of reaction species is 
proportional to the enthalpy gradient; that is, 


q 


c 

z 


q 


c 

r 


dh 

dz 

dh 

dr 


From conservation of energy across the near normal portion of a strong spherical shock 


Thus, 






c 2 H R s 
q r ^ r 


By comparing terms that appear in the energy equation, it is found that 



8qC/ 9z 



« 1 
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Thus, the terms containing q£ can be neglected in the formation of the stagnation-flow 
model. 

If the shock layer is optically thick, that is, the photon mean free path is very 
small compared with the shock standoff distance, the radiation flux terms take on the 
same form as the conduction term (according to the Rosseland or diffusion approximation) 
and 8q^/8r + qR/r may be neglected. On the other hand, if the gas is not optically 
thick, this simple order of magnitude analysis no longer suffices because the divergence 
of the radiant energy flux depends not only on local conditions, but on conditions through- 
out all the shock-heated gas which can be seen by an observer located at the point in 
question. 

Calculations were made of the divergence of the radiant flux for a gray isothermal 
gas in a shock layer formed by two concentric spherical surfaces with a standoff distance 
to shock radius ratio of 0.05. A sketch showing the volume of gas which contributes to 

the radiant flux at a point on the 
stagnation streamline is shown in 
sketch (d). 

The largest difference between 
this result and the divergence of the 
radiant flux for a plane -par all el layer 
occurred adjacent to the wall for an 
optical thickness of about 0.13. The 
difference amounted to 2.2 percent 
of the value for the plane-parallel 
layer. 

A second set of calculations 
was made to determine the effect of 
a nonuniform enthalpy distribution, 
in the lateral direction, on the mag- 
nitude of the divergence of the radi- 
ant flux. The enthalpy distribution 
was given by 



h(r) 



1 ° 


(r £ \/2Rs) 
(r > \J2R S ) 


(34) 


This expression approximately corresponds to the enthalpy distribution in the shock layer 
about a spherical body. The absorption coefficient was assumed to vary as the third power 
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of the enthalpy (this variation is consistent with the correlations of the optical properties 
of air to be discussed in a subsequent section) and the shock standoff distance to shock 
radius ratio A/r s was chosen to be 0.05. A comparison of calculations for a plane- 
parallel layer in which the enthalpy was assumed to vary according to equation (34) and 
of calculations for a plane -parallel layer in which the lateral enthalpy distribution was 
uniform (that is, h(r) = h(0)) indicated that the largest difference in the magnitude of the 
divergence of the radiant flux occurred for a shock-layer optical thickness of about 0.1. 
This difference amounted to 2.8 percent of the value for the uniform distribution. These 
results are independent of the value of h(0) . 


Since the errors in the divergence of the radiation flux for the one -dimensional 
shock layer due to the separate effects of geometry and nonuniform lateral distribution 
of enthalpy are small, their combined effect should be given approximately by the sum of 
the separate effects. That is, the maximum error due to the combined effects of geom- 
etry and nonuniform lateral distribution is probably not much greater than 5 percent for 
A/R s = 0.05. This statement, of course, does not imply that the final results for the 
enthalpy (for example) would be in error by 5 percent but only that one term in the energy 
equation is in error by 5 percent. In any event, the results of the calculations mentioned 
are considered to give sufficient justification for choosing the plane -par all el layer as a 
model for the stagnation region of a blunt body. ^ 


The expression for the radiant energy flux is more seriously affected by the plane- 
parallel layer assumption than is the divergence. For example, Koh (ref. 21) has com- 
puted the radiant flux at the wall for an isothermal shock layer formed by two concentric 
spheres. For a shock standoff distance to body nose radius ratio of 0.05 and a vanishingly 
small value of optical thickness, the result is about 17 percent less than that for a plane- 
parallel isothermal layer of equal optical thickness. This difference decreases with 
increasing optical thickness. Koh also computed the effect of nonuniform lateral enthalpy 
distribution by using an assumed enthalpy distribution similar to that given by equa- 
tion (34). He found that the flux at the wall for the nonuniform distribution was about 
1.5 percent less than that for an isothermal layer for a shock standoff distance to shock 
radius ratio A/R s = 0.05 and a vanishingly small optical thickness. As expected, the 
difference decreases as the optical thickness increases. 


It is apparent from Koh’s results, that an accurate estimate to the rate of radiant 
heat transfer to the stagnation point cannot be obtained through the use of the plane- 
parallel layer approximation unless some correction factor, which takes into account the 
actual geometry of the shock layer, is employed. However, because this investigation is 
concerned with obtaining a general understanding of the problem of radiating shock layers 
rather than with specific numerical results, such a correction factor will not be used / 
herein. 
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At this point, it is convenient to introduce the variable transformation 



( 35 ) 


The new variable 7) is often called the Dorodnitsyn variable. Under this transformation, 
the normal and tangential velocity components become 

w = -2f (if) — (36a) 

dri 

u = rP(rj) (36b) 


The two momentum equations (eqs. (16) and (17)) take the form 

(pf”(^’ + 2i( V ) f»fo) - [f*(^ 2 = 


(37) 


and 


- 2 p 


p( f '(?7) - 1 p'(v) f(t7)) + p (p' - 1 p) ^ p’(? 7) f(r?) 


2i(v) i'(v) 


+ | p'Mpfaj] 2 = +| PP’fa) 


(38) 


An order of magnitude analysis of equation (38) indicates that p '(rf) is order x 
or N^, whichever is larger. Since both y and are very small compared with 

unity, equation (38) will be replaced by the simple approximate expression 

p’O?) = 0 (39) 

Thus, the pressure is a function of r only. In particular, the strong shock relations for 
the near-normal portion of a spherical shock give 


P( r ) = PooW^ 1 " x) 


(£) 2 


o(r 4 ) 


To first-order in r 


- iE. = -2 p w 2 (l - y)— 

r 9r ^oo TT ooV a; 9 


Re 


so that equation (37) becomes 


[ppf"(7?)] + 2f(?7) f ”(??) -[f’(?7)] 2 = ~ — ' ■ °° g — 


pRi 


(40) 


(41) 
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Under the foregoing assumptions and the coordinate transformation, the energy 
equation (18) becomes 

-2f(rj) h’fa) +q’( 77 ) = 0 (42) 

The boundary conditions are 


f( 0 ) = 0 
f'( 0 ) = 0 




f 


( 7 ? a) = 


Wo 

R c 


h( 0 ) = h w 


h w = \ 


w 2 

2 00 


(43) 

(44) 

(45) 

(46) 

(47) 

(48) 


where 

P A 

V A = J o P dz (49) 

The heat flux term q( 7 j) in the energy equation (42) is composed of a combined 
conduction and diffusion term 


Q C 0?) = -pkT* ( 77 ) = h’ (??) 

Np r 


and a radiation term 


(50) 


q R (? 7 ) = -277 






B X (tx) E 2 (tx) dt A 


dX 




(51) 


This radiation term is representative of the case of a plane-parallel geometry with a 
transparent wall (shock) and a cold wall, which reflects diffusely and independently of 
wavelength a fraction r w of the incident radiation, separated by a nonscattering nongray 
gas. The variable t x is called the "optical path length" and is defined by the expression 
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nZ nrj 

T X - J 0 *>«x dz = J 0 K x d ” 


( 52 ) 


Expression (51) was specialized from the more general expression of Goulard (ref. 22). 
Goulard derived the expression for the radiant flux in a plane-parallel geometry with arbi- 
trary reflecting, absorbing, and emitting walls separated by a nonscattering, nongray gas. 
His expression was restricted to isotropically emitting and diffusely reflecting surfaces 
and hence, so is equation (51). 

The first term in equation (51) represents the radiant energy flux passing through 
the plane = Constant and which originated in the region between this plane and the 
shock at This radiant flux has been attenuated by partial absorption in the 

intervening gas. The second term represents the radiant energy flux passing through the 
plane = Constant and which originated in the region between this plane and the wall 
at = 0. This flux has also been attenuated by partial absorption in the intervening gas. 
The last term represents the radiant flux passing through the plane t x = Constant and 
which was reflected from the wall and attenuated by the intervening absorbing gas. 


Substituting the expressions for the energy flux equations (50) and (51) into the 
energy equation (42) gives 


2f ( 77 ) h f ( 77 ) 


h' (77) 

Np r 


•+ 


ip] = 0 


(53) 


where the divergence of the radiant flux is represented by the integral term 



= -47TKp(77) 13(77) + 2?r 



K k W) ('//') e/ 


t aW 



d77’ 


+ 2r w E2 



K x ( 7? f ) B x ( 77 t ) E 2 (t a ( 77')) d77' 


dX 


(54) 


The final step in the derivation is to reduce the equation to nondimensional form. 
For this purpose, the following set of nondimensional quantities is introduced 
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V = P S A pj}, *(V) = | PooWjtf), h(7?) = i W 2 h(rj) 




pm _ (W 
N Pr \ N Pr,s y 

,4 


P^h), p/i = (p s M s p 2 00, P -1 = pg^gte) 


1 &] = 2oT s K P,s , K \W = /c PjS K x (i 7), r x { 77) = k p r x (r 7) 


L 


dm 

\ \hJQ. Wl 

\C.^> : ^5 0 (A -(Wavx OOv 


(Equations continued on next page) 


J 





^S C D S Poo W oo A A l 

^Pr,s ^ 1 ’ ^Re,s JT ~ > -^Pe,s - ^Pr,s ^Re,s > 


P U- a> A A 


£ 


/tA- 




!)^ 


n Bo = 


4oT 


4’ k P “ P s K P,s A A’ 


e = 


N 


Bo 


(55) 


4 <H, 




VJ 


The subscripts °° and s indicate conditions evaluated in the undisturbed free stream 
and immediately behind the shock, respectively. The quantity A ^ is the shock- standoff 
distance for the nonradiating (or adiabatic) shock layer. The property variations repre- 
sented by *2^, jf 2 > are functions of h only as the gas is in local thermody- 

namic and chemical equilibrium and the pressure has been assumed constant throughout 
the stagnation region. The quantities Np r s , N Re s , and Np e s are the Prandtl, 
Reynolds, and Peclet numbers, respectively, based on conditions immediately behind the 
shock. The parameters Ng 0 , kp, and e are the Boltzmann number, the Bouguer num- 
ber, and the radiation cooling parameter, respectively. These parameters are fundamen- 
tal to the study of radiation gas dynamics and have been discussed by a number of inves- 
tigators. (See, for example, refs. 23 and 24.) 

Substituting these nondimensional quantities into equations (41), (43) to (48), (53), 
and (54) yields the nondimensional system governing the flow in the stagnation region of 
a blunt body traveling at hyper sonic speeds 


f (77) h' (77) 




(56) 


2N, 


e,sL 


~ ~ f 2 

7 2 (h) f"( 5 j) + 2f (77) f (77) - [f (77)] + a 2 3 3 (h) = 0 


f (0) = 0 


f (0) = 0 



1 



a 

\/2x(l - X) 


(57) 

(58) 

(59) 

(60) 

(61) 


19 



(62) 


h(0)=h w = r 2L 

h(^ A ) = 1 


(63) 


where rj^ is the value of the nondimensional Dorodnitsyn variable at the shock. This 
quantity is determined from the expression 


p^A_ 

= h(rf) drj 
Jn 


_A 

a A ^0 


(64) 


The integral term I £ 77 ] is given by the expression 


l[v] = -2 k p (v) B(~n) + kp J K x (r]) I J K^v) Ejfk p - T x (ri) 


d t? 


2r w E 2 (kpT A ()7)) k x (v) B x (v) E 2 (kpr^( t)')) dr? 


dX 


(65) 


Later in this paper, it will be convenient to express the energy equation in terms of 
the optical path length as an independent variable. In both cases the optical properties of 
the gas are assumed to be independent of wavelength. In this event, the energy equation 
(less the thermal conductivity term) becomes 

f (r) h (t) + el[rj = 0 (66) 

where 



+ 2r w E 2V k P T) 


i( k p f ) I 


A B(t) Er 


(kpt) 



(67) 


Throughout the remainder of this paper, the bars over the nondimensional variables 
are dropped. This procedure should not lead to any confusion because only the nondimen- 
sional form of the governing equations are employed. 


Divergence of Radiant Flux 

The nondimensional form of the divergence of the radiant flux is 
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( 68 ) 


eI H = -2e/c p (?7) B(v) + ek p 



K x (r]') B x (t 7') Ej^kpj 7^(77) ~ \(r]') |) d77 f 


+ 2r w E 



k. P 

r 

K \(V') B x (t7') E 2 (k p t x (? 7’)) dry’ 


dX 


The first term on the right-hand side of this expression is the local emission term which 
represents the rate at which energy is emitted per unit volume of gas at the location 77. 
The integration over all wavelengths X has been performed for this term with the aid of 
the definition of the Planck mean-mass absorption coefficient. The second and third 
terms represent the rate at which radiant energy is absorbed per unit volume at the 
location 77. 

It is the presence of the second and third terms which so greatly complicate the 
radiation problem. These terms are integral expressions. In addition, their presence 
makes it impossible to define a wavelength averaged absorption coefficient by which the 
wavelength dependence might be eliminated. The importance of these terms is indicated 
by the magnitude of the Bouguer number kp which is the ratio of the shock standoff 
distance for a nonradiating shock layer to the photon mean free path evaluated at condi- 
tions immediately behind the shock. 

Tjie ^adia tionj; ooling parameter e is a ratio of the rate of energy loss per unit 
area by radiation from both sides of a nonabsorbing isothermal layer of gas of thick- 
ness A a to the rate at which kinetic energy enters the shock layer per unit area of 
shock surface. Alternatively, the parameter e may be interpreted as the ratio of the 
radiationless standoff distance to the decay length where the decay length is the length 
required by an element of gas to lose all the energy it possessed upon emerging from the 
normal shock if it loses this energy by radiating (without reabsorbing) at a constant rate. 
This parameter modifies the entire radiation term and thus, acts as a measure of the 
relative efficiency of radiation compared with convection as energy transport mechanisms 
within th e shock layer. In addition, the surface reflectivity r w and the wavelength 
dependence of the absorption coefficient influence the character of the radiation terms and 
are considered as parameters in this study. 

Most investigators who have studied problems in which a term similar to 1Q7J 
appears have assumed that the gas and its surroundings are gray; that is, the optical 
properties are independent of wavelength. This assumption allows the integration over 
frequency to be performed analytically. Accurate results can be achieved in the two 
extreme cases of optically thin « l) and optically thick (t^ » 1) gases. When the 
gas is optically thin at all wavel engt hs , the gray absorption coefficient is correctly given, 
by the Planck mean mass absorption coefficient 
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dA j VlcvV^c^* 


(69) 
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= j" 


where is the monochromatic mass absorption coefficient and the weighting func- 
tion Bx is the Planck black-body function. When the gas is optically thick at all wave- 
lengths, the gray absorption coefficient, in the interior of the gas , is correctly given by 
the Rosseland mean mass absorption coefficient 


K 


R ' 




~1_^X 
X 9X 


dX 


(70) 


Near a rad iation bound ary or in regions of rapid (w ith resp ect to the optical path length) 
variat ions in thermodynamic properties the Rosseland mean is not valid. At intermediate 
values of optical depth, no single mean absorption coefficient, which depends only on local 
thermodynamic conditions can be defined. In fact, as has been pointed out by Krook 
(ref. 25), it would be necessary to define an infinite number of such mean coefficients. 

This restriction, of course, does not preclude the possibility of defining approximate mean 
coefficients under these conditions. 

Stone (ref. 26) introduced a model in which the monochromatic absorption coeffi- 
cient was a step function of frequency, the size of the steps being independent of the geom- 
etry or thermodynamics of the system. By means of this method, the integral over all 
wavelengths is reduced to a finite series. Carrier and Avrett (ref. 27) considered an 
absorption coefficient with only two steps, one of which was very much larger than the 
other. Both of these papers were concerned with Milne's problem of a stellar atmosphere 
in radiative equilibrium. Lick (ref. 11) and later Greif (ref. 12) studied the problem of 
one-dimensional energy transfer between two walls separated by a radiating and con- 
ducting gas. A picket fence model, which is a specialization of the step function model, 
for the absorption coefficient was used. Krook (ref. 28) derived expressions by means 
of the Poincard-Lighthill-Kuo perturbation procedure for a slightly nongray gas. The 
solution represents a perturbation to the gray gas solution. Rhyming (ref. 29) considered 
wave propagation in a simple dissociating flow of a radiating gas where the absorption 
coefficient was given as a Gaussian function of the frequency. 

However, even with these simple models for the absorption coefficient, the term 
l[Vj retains an integral character and the solution to the set of equations is still very 
difficult to obtain. Numerical procedures are extremely tedious. As a result of this 
difficulty, several approximate analytical methods have been derived in order to reduce 
this term to purely differential form. One such technique is the Milne-Eddington approxi- 
mation (ref. 30), the derivation of which has been based on physical considerations, but 
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which may also be thought of as a substitute kernel approximation (ref. 31). The integral 
terms can then be eliminated by means of a double differentiation (for a gray gas only). 

| Of course, this increases the order of the differential equation by two. This technique has 
| been used by a number of authors in the study of the dynamics of radiating gases. (See, 

1 A 

I for example, refs. 32 and 33). Barbier (ref. 30) introduced the method of expanding the 
| source function in a Taylor series about the zero of the argument of the exponential 
1 integral kernel, E^| T^(r}) - t^(?? t ) j Because the kernel function has a logarithmic 

' singularity at the zero of its argument, the integral over the first few terms of the series 

8 

Sshould provide a good approximation. The resulting integrals can then be evaluated ana- 

if 

|lytically and the equation becomes purely differential in character. Yoshikawa and 
|Chapman (ref. 8), Thomas (ref. 35), and Viskanta (ref. 36) all used the method of Barbier 
|to different degrees of approximation. When more than the constant term in the Taylor 
^series is retained, it may become necessary to introduce additional boundary conditions. 


'In fact, it may not be enough merely to specify a new condition, it may be necessary to 

§ / 3 

| modify the existing conditions as well so that in the limit as the parameter N = k/Cp^oTg 

I (for example) tends to zero, the solution will approach the proper pure radiation solution. 

f 

E This parameter, which appears in the literature concerning energy transport by radiation 
and conduction, represents the relative importance of conduction compared with radiation. 
i-When N tends to zero, radiation is the dominant mode of energy transport. 


The diffusion approximation for optically thick gases has been used extensively in 
| astrophysics and gas dynamics. Probstein (ref. 37) has shown how to extend the useful- 
\ ness of this approximation to gas layers of finite optical thickness by means of radiation 
j slip-boundary conditions. It is not at all clear, however, that these slip conditions can be 
used in the problem of this paper because of the presence of the convection term 
i pw dh/dz. 

i 

j The optically thin approximation of hot gases, in which absorption is neglected in 

\ comparison with emission has also been used extensively in gas dynamics. As Thomas 

I (ref. 35) has pointed out, this approximation is not valid in those portions of the gas which 

J 

sare considerably cooler than the remainder of the gas. 

/ 

| In this paper, the integral term l\jf\ is reduced to algebraic or differential form 
fthrough the use of various approximations similar to those described. The manner in 


I 


| ^Traugott (ref. 34) has introduced a "method of moments” in order to reduce the 
^integral term to differential form. This method may be taken to any degree of approxi- 
mation desired (not without a considerable sacrifice in simplicity, however). The first 
^approximation is identical to the Milne-Eddington approximation. Traugott's higher 
{approximations can also be obtained by a substitute kernel method similar to that of 
[Krook (ref. 31 ). 
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which this reduction is to be accomplished will depend on the order of magnitude of the 
parameters e and kp and is discussed in detail in the subsequent sections. 

The Inviscid Shock Layer 

As was pointed out in the Introduction, the studies of this paper will be concerned 
only with those cases for which the thicknesses of the wall boundary layers due to the 
presence of viscosity and thermal conductivity are very much less than the shock stand- 
off distance. For a nonradiating gas, the shock layer can be separated into an outer 
inviscid and nonheat-conducting region and an inner viscous and heat-conducting region 
or boundary layer. Considerable simplification will result if a similar separation can be 
achieved in the case of a radiating gas. As will be shown, such a separation can be 
obtained when the boundary is either optically thin or optically thick. Only the former 
situation is considered herein. The method of separation follows the procedures delin- 
eated by Van Dyke (ref. 38). Mathematical details are presented in appendix A. 

It is shown in the appendix that the significant parameter which determines the 
extent of the boundary layer is the inverse square root of the Peclet number, (Np e )"'^. 
The zero-order in (Npg) - *^ system of equations which governs the flow in the inviscid 
region is 


foM h o(^) + eI oH = 0 

(71) 

2f 0 (r?) fg(7?) - [fg(i])] 2 + a 2 3 3 (h) = 0 

(72) 

o 

II 

3 

'"o 

(73) 

f o0k) = 1 

(74) 

f „ \ 2 ^A a 

°™ XR S ^(ir- x ) 

(75) 

h o{ v *) = 1 

(76) 


The dependent variables f Q (?j) and h Q (? 7 ) are the asymptotic values of the non- 

dimensional stream function, and h(? 7 ), the nondimensional enthalpy, respectively, as 
fNpe)~ 1/72 approaches zero. 
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In the boundary layer, the zero-order system of equations is 


i 

[(?l(i 0 ) tyO] + g 0 (S) fy*) + eJ oH = 0 

(77) 

I 

si 

f 

1 

2N Pr,s[7 2 ( i o) g 0 ( ^] + 2 g 0 (?) g 0 ( ^ " [ g 0 ( ^ + = ° 

(78) 

i 

tf 

g 0 (0) = 0 

(79) 

J 

! 

V 

gq(0) = 0 

(80) 

i, 

1 

lim g Q (0 = fUO) 

(81) 
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£-»oo J 


- 

s 

!' 

i 

ig(0) = h w 

(82) 

- 

{ 

lim i Q (|) = h Q (0) 

(83) 


The independent variable £ is the "stretched" boundary-layer coordinate defined by the 
relation 

{=( N Pe) 1/Z 1 < 84) 

The dependent variables g(£) and i(£) are defined by the expressions 

i(l) = Hv) 
g'U) = i'(v) 


(85) 

( 86 ) 


■ 1/2 


3; in the boundary layer as (Np e ^ “ approaches zero. 

I; The quantities I q[V] and J 0 (jQ are terms of zero-order in the expansion of 

I [ 773 , the divergence of the radiant flux vector, for the inviscid and boundary-layer regions, 
respectively. The derivation of these quantities is presented in appendix A. The results 
are: 
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IqW = -2/c p [h 0 (^)] B[h 0 (77)] 


and 



J o[d “ - 2 -p[Vd B (V^j 

E 1 ^k p r x (77')) d77* 


dr/'J dA (88) 

Th e int egrals which appear in the second of these expressions are definite integrals. 
Consequently, the system of equations governing the flow in the boundary layer is a purely 
differential system. 

It must be realized that expressions (87) and (88) are restricted to the case of an 
optically thin boundary layer. It is only in this case, and the case for which the boundary 
layer is optically thick, that a complete separation between the inviscid region and the 
boundary layer can be achieved. At intermediate values of optical depth, the integral 
term Iq^V] i- s a function of the enthalpy distribution in the boundary layer in addition to 
being a function of the enthalpy distribution in the inviscid region so that the equations in 
the inviscid region and the boundary layer are coupled. The i nfluence on the inviscid 
region of radiation from an optically thick boundary layer cannot be neglected. However, 
most of this radiation originates at the outer edge of the boundary layer. The boundary- 
layer solution in this region is constrained by matching conditions to approach asymptoti- 
cally the value of the inviscid solution at the wall. Hence, the radiation contribution to 
the inviscid region from the boundary layer can be obtained from the inviscid solution at 
the wall, and the inviscid solution remains uncoupled from the boundary-layer solution. 


k P J 0 Vo A ^[ h ° (7?,) ] B *[ h o (77 ' } ] 


2r w A K A[ h 0^[] B a 1VV] E 2(kpT A (T7’)) 
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This restriction to an optically thin boundary layer is not so severe as it might 
first appear because the Planck mean optical thickness of a boundary layer in whic h th e 
absorption coefficient is the same order of magnitude as it is for shock-heat ed air will 
not exceed about 0.1 at any altitude (above 20 km) and velocity (up to 15 km/sec) for a 
shock radius of 1 meter or less. In fact, the Planck mean optical thickness of the bound- 
ary layer will be less than 0.1 at that point of the trajectory of a Martian or lunar return 
vehicle with a shock radius of about 1 meter for which heating is a maximum even if the 
absorption coefficient in the boundary layer is 100 times that of shock-heated air. That 
the optical thickness should be small is not so difficult to see when it is realized that both 
the optical-path length and the boundary-layer thickness decrease rapidly with decreasing 
altitude. Thus, at low altitudes where the optical-path length is small and the shock layer 
may be optically thick, the boundary-layer thickness is very small. For larger objects, 
the boundary layer need not be optically thin at the lower altitudes because the boundary- 
layer thickness depends on the size of the object whereas the optical-path length does not. 

These conclusions regarding the optical thickness of the boundary layer generally 
concur with the observations of Fay, Moffatt, and Probstein (ref. 10). Henceforth, the 
discussions of this paper will be limited to the case of an optically thin boundary layer, 
and radiation from this boundary layer will be considered to have no effect on the solution 
in the inviscid region of the shock layer. 

If the inviscid system of equations (eqs. (71) to (76)) is solved for the nonradiating 
case (e = 0) along with condition of equation (64), the ratio of the shock standoff distance 
to the shock radius is given by the expression 


a A 

R« 


(89) 


1 + )/2x(l - X) 

Hayes (ref. 39) obtained the same result when the shock and body surfaces near the stag- 
nation point are concentric spheres. When the shock and body surfaces are not concentric 
(that is, R s ± + A), condition (89) is still approximately true over a wide range of 

body shapes. (See, for example, refs. 40 and 41.) With this result 


a = 


2^ x (l-x) 


1 + \/2x(l " X) 

This value for a, the constant appearing in the momentum equation (72), will be used 
throughout the remainder of this investigation. 

Thermodynamic and Optical Property Correlations 

In order to achieve meaningful results, an attempt has been made in this paper to 
use simple yet physically reasonable approximations to the thermodynamic and optical 


(90) 
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properties of high temperature gases. In particular, correlation formulas were derived 
from the existing store of information about equilibrium air. The thermodynamic proper- 
ties were obtained from reference 42, for temperature up to 100 000° K and pressures 
from 10~® to 10^ times atmospheric. The optical properties were obtained from a variety 
of sources which will be cited later. 


It was noted from the data of reference 42, that both the density and temperature 
could be approximately represented by functions separable in the variables pressure 
and enthalpy; more specifically in the form (^p/Pg^ n f(h/RTg^. Plots of the functions 
f(h/RTg) for the density and temperature at various pressure levels are presented in 
figures 1 and 2, respectively. It is apparent from these plots that the density and tem- 
perature can adequately be represented by the expressions 
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A number of investigators (see, for example, refs. 43 to 47) have calculated the 
radiant properties of equilibrium air for temperatures up to 25 000° K and for densities 
from 10"® to 10* Amagats. Because of the extremely complex nature of these calcula- 
tions, the many physical processes which produce radiation, and the uncertain knowledge 
of cross sections and transition probabilities the scatter among the various calculations 
is often quite large. Some of the results for the Planck mean mass absorption coeffi- 
cients are presented in figure 3. 

A correlation formula can be obtained from figure 3 by approximating the curves of 
the variation of log^Q p/Cp with logjg T with straight lines. The resulting formula is 



This formula is not a particularly convenient form for use in the calculations of this paper. 
It is more desirable to express the Planck mean mass absorption coefficient Kp in terms 
of the pressure and enthalpy. This expression was determined by crossplotting the loga- 
rithm of the absorption coefficient data shown in figure 3 against the logarithm of the tem- 
perature at constant pressure. Straight lines were then fitted to the resulting curves. 
Finally, the correlation formulas (eqs. (91) and (92)) for density and temperature were 
used to obtain the formula: 
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K-p = (3.82 X 10~ 9 )(§ 


+1.425-0.44 lo gl 0 (^J ^ 3.55-0.24 log 10 (L 


This formula is valid for temperatures up to 20 000° K at the higher pressures 
(p/Pq = lO - "*" to 10 1 ) and to somewhat lower temperatures at the lower pressures (for 
example, when p/Pg = the maximum temperature at which the formula is valid 

is 15 000° K). 

At temperatures above about 20 000° K, the information about the radiative proper- 
ties is not so comprehensive. Most of what exists consists of Planck and/or Ross eland 
mean absorption coefficients for continuum radiation. Line radiation is neglected. At 
these high temperatures, the radiation consists of spectral lines of the various ions which 
may be appreciably "Stark broadened" at high electron densities, and a continuum due to 
free bound and free-free transitions of electrons in collisions with the ions. S ince jthe 
integrated line emission is proportional to the ion density whereas the continuum emis- 
sion is proportional to the product of ion and electron densities, the ratio of the latter to 
the former increases with increasing density. Thus, at the higher density levels, the 
continuum calculations may be adequate. 

The results of calculations of the Planck mean volume absorption coefficient for 
temperatures from 23 200° K to 100 000° K are presented in figure 3. These results 
were obtained from the paper by Armstrong et al. (ref. 45). A correlation formula was 
obtained from these results crossplotted in terms of constant density. The resulting 
formula was 


/c p = (2.95Xl0 19 )^' 


-0.055+0.011 logjg^-j -1.217+0.065 log 10 ^ 


, cm^/gm (95) 


This formula can be used in the range of temperatures from 30 000° K to 100 000° K. 


THE SMALL -PERTURBATION SOLUTION ^ 


The Conventional Method 


! There is a flight regime of considerable importance in which the radiation cooling 

parameter e is very much less than unity. In this regime, the energy transferred by 

^The material in the corresponding section of the thesis has more recently appeared 
!", in reference 47 along with a considerable number of calculations for an improved nongray 
air absorption coefficient model. This report will present a summary of what is covered 
! in the thesis and the reader is referred to reference 47 for details. 
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radiation is small compared with the influx of kinetic energy across the bow shock, and it 
would be reasonable to expect the flow properties to be only slightly perturbed from the 
radiationless case. Lunev and Murzinov (ref. 4) and Goulard (ref. 5) took advantage of 
this property and developed what amounted to first-order perturbation solutions of the 
temperature distribution in the inviscid region of a transparent, gray gas layer. In theses 
papers, simplifying assumptions concerning the gas properties and flow model have been 
included. 

In this section, the perturbation solutions are generalized to include nongray gases 
with arbitrary thermodynamic and optical properties. These solutions are not limited to 
shock layers of small optical thickness. Also, the solutions are extended to second order. 
As will be shown, the second-order solutions can be very important when the absorption 
coefficient varies rapidly with temperature. 

The integro-differential system which governs the flow in the inviscid region of the 
shock layer is 


f(v) h '(v) + el [77] = 0 

(96) 

2f(rj) f"(v) ~ [f '(V)] 2 + a 2 h(7?) = 0 

(97) 

f(0) = 0 

(98) 

f K) = 1 

(99) 

h ( ? ?A) =1 

(100) 

fin \ _ 2 \ _ __ a 

v A ) x\R s j 

(101) 


Here f ( 77 ) and h (rj) are the nondimensional stream function and enthalpy, respectively. 
The quantity t? a is the value of the Dorodnitsyn coordinate at the location of the shock. 
The constant a can be expressed in terms of x> the density ratio across the shock, 
through expression (90). When the radiation cooling parameter e is very small, the 
integral term in equation (96) becomes of only secondary significance throughout most of 
the domain of the problem with the obvious exception of the region 77 » 0 where f(rj) ~ 0. 
(The difficulties presented by this exception are discussed later.) Neglecting the integral 
term I^ryJ reduces the problem to one in which radiation does not play a part. If, as 
expected when e is small, the presence of radiation only slightly influences the solution 
one can, to reasonable accuracy, evaluate l(VJ by using the radiationless solution for h 
so that equation (96) becomes purely differential. Thus, when the small -perturbation 
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procedure (which roughly proceeds in the manner outlined above) is applied to this prob- 
lem, the integro-differential system is simplified to a purely differential system. In 
addition, as a result of the nature of the lowest order solution for the enthalpy distribu- 
tion, the two differential equations become uncoupled and can be solved independently. 
Hence, it becomes possible to obtain analytic solutions to any order of approximation to 
the flow in the inviscid region of the shock layer. Details of the derivation of these solu- 
tions are presented in reference 47. 

The zero-order, or radiationless, solution is simply 

h 0 (??) = 1 (102) 


i Q (v) = (1 - a)?? 2 + a?7 


(103) 


The first-order solution, which represents the effect of radiation when it is assumed that 
the emissive power of the gas is independent of temperature, is 

I 0 [x] dx 


hi (V) = J 


?7 (1 - a)x 2 + ax 


(104) 


flM = “ a) 7 ? + a) ^ l(X) -- 2 -dx - V 2 j* ~ a | X - + ~ - ^(x) dx\ (105) 

2 V- JJ 0 [(l-a)x + a] 2 J ?7 x2[(l - a)x + a| 2 J 

Here x is a dummy variable of integration. The quantity Io[V] is given by the formula 

r oO] = "J 0 + i 1 - r 0 ,x) E 2[ k X^} ^ ( 106 ) 

The notation has been simplified somewhat in this expression by omitting the argument hq 
in the terms k x and and by introducing the quantities 


k x = kpK x 


Also 


r 0,X = r w[! - 2E 3( k x)] 


$l(v) = - \ a 2 h 1 (?j) 


(107) 

(108) 

(109) 


The second-order solution takes into account the change in gas properties with 
changes in enthalpy. This solution is 
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( 110 ) 


f 2 W = -(1 - a) 2 ^ ^ - |j2(l - a)ij + a] ^ — 


V $o(x) 


dx 


J(1 - a)x + aj^ 


+ h 2 L 


2(1 - a)x + a 
x 2 j(l - 2)x + aj 


^ $ 2 (x) dx 


(HI) 


where 


IiH = h l^ J 0 (^ B x{ E 2[V 1 - + (! - r O,x) E 2[V 


- 2k x B a dX 


I 0 k x( K x B x j 0 h iW E i[ k xh " X [| ** + "x^^ifxd - ^ ^ h lW * 
(! - r 0,x) E l[ k X ?7 ] JJ 7 h l (x) + r l,X E 2(V]j dX 


( 112 ) 


Here the argument hq is omitted in the terms and and the quantity r ^ ^ is 

defined by the expression 


r l,X= 2r w 


• r 1 r 1 

l B xJ„ h X x ) E 2( k x x )*‘ + ' ; x B x E 2( 1[ x)J 0 h l( x) dx + , A1 « x B x E 2 (k x ) 


(113) 


Also 


$ 2 0?) = + 1 f i(^?)J 2 - | a 2 h 2 (?7) 

The quantities q> 1’ and 2 are &i ven by the formulas 


1a,0 = 1 

1 pi ^l(x) 

1 A < = M ffdx 

A)1 2 J 0 \ii - a )x + a] 2 


(114) 

(115) 

(116) 
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(117) 


’>A,2 = < 1 - a >4,l + f.f„ 1 


$?(x) 


dx 


0 [(1 - a)x + aj 


It can be seen upon inspection of relation (112) that a large value of the rate of 
change of the Planck mean absorption coefficient with enthalpy will lead to large values 
of Thus, at shock temperatures of less than about 30 000° K, for which the 

absorption coefficient does vary rapidly with enthalpy, the second-order solutions can 
become more important to the overall solution than their order in e might at first 
indicate. 


The Poincare-Lighthill-Kuo Method 

Inspection of the expressions (104) and (110) shows that the first-order solution for 
the enthalpy distribution has a logarithmic singularity at the point 77 = 0 and the second- 
order solution has a singularity of greater strength at this point. As a consequence, the 
assumed expansion diverges as the origin is approached and the small perturbation solu- 
tion is not uniformly valid throughout the domain of the problem. This divergence can 
lead to serious errors in the calculation of the radiant heat flux to the wall because those 
regions close to the wall, in which the largest errors occur, are given the most weight in 
the calculation. This condition is particularly true for shock layers which are not opti- 
cally thin. Additional difficulties are encountered when an attempt is made to specify the 
proper outer boundary conditions for the viscous boundary-layer equations. In classical 
boundary-layer theory, the outer boundary conditions are obtained from the values of the 
outer (or inviscid) solution at the wall (77 = 0 in this problem) . Because of the divergence 
of the outer solution, no finite value exists at 77 = 0. 

In this section, the Poincare-Lighthill-Kuo (P-L-K) perturbation of coordinate pro- 
cedure^ (ref. 48) is used to obtain a solution which is uniformly valid over the domain of 
the problem. The details of the application of this method to the problem of this paper 
are presented in reference 47. This method utilizes a coordinate transformation in the 
form of a perturbation expansion of the coordinate to remove the singularity (which caused 
the divergence of the conventional solution) from 77 = 0 to a small negative value of 77 
which lies outside the domain of the problem. The P-L-K expansions are 

77 = X + 6 77 1 (x) + . . . (118) 

h(77;e) = 1 t.q(x) + eh^(x) + . . . (119) 

f(??;e) = f q( x ) + ef iM + • • • (12°) 

D Variously called the P-L-K method, the P-L method, Lighthill's technique, the 
method of strained coordinates, and the method of perturbation of coordinates. 
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where x is the coordinate in the transformed plane, and the asterisk superscript has 
been used to differentiate between the coefficients in the P-L-K expansion and the coeffi- 
cients in the conventional expansion. Pritulo (ref. 49) has derived a general relation 
between the P-L-K and conventional coefficients. Adapted to this problem, the relation- 
ships become 


h Q (x) = M x ) 

(121) 

hj(x) = hj(x) 

(122) 

fg( x ) = f o( x ) 

(123) 

f*(x) = f^x) + 77* (x) f q(x) 

(124) 

*, , " h 2 W 

(125) 


The second-order term h^x) introduces the effects of variable thermodynamic and 
optical properties; thus it is apparent that these effects are contained in the first-order 
P-L-K solution. 

A comparison of the P-L-K and conventional perturbation solutions for the enthalpy 
distribution for a constant-density transparent shock layer is presented in figure 4. The 
divergent character of the conventional solutions is apparent. Also shown in this figure 
is the exact analytic solution which can be obtained in this simple case. The formula for 
this exact solution is 

1 

y-l 

> (126) 

where y (the exponent in the correlation formula /CpB = h^) was taken to be 6 and the 
constant a (which appears in the momentum eq. (97)) was taken to be 0.5. The good 
agreement between the P-L-K solution and the exact solution indicates that the accuracy 
of the P-L-K solution is probably second order in the radiation cooling parameter e 
throughout the domain except in the immediate neighborhood of the wall. It is clear that 
quantities such as the radiant heat flux at the wall, which depend upon an integration over 
the enthalpy distribution, will be considerably more accurate if the P-L-K solution rather 
than the conventional perturbation solution is used. 

It should be noticed that the P-L-K solution does not lead to zero enthalpy at the 
wall as the exact transparent solution does. The reason for this disparity can be found 
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in the fact that the coordinate stretching displaces the boundary with regard to both the 
energy and momentum equations but not by a uniform amount. Thus, a physical inter- 
pretation of the first-order P-L-K solution is that the normal velocity of the flow at the 
boundary for the energy equation is not exactly zero, and a particle approaching this 
boundary will reach it in a finite time before losing all its energy by radiation. 

It can be shown that since the expected error in the Dorodnitsyn coordinate 77 in 
terms of the stretched coordinate x is order and since the gradients in h^x) 
are very large in the vicinity of the wall, the difference between the P-L-K and exact 
solutions at the wall lies within expected limits. Convergence to the correct solution may 
be attained with the addition of higher order terms to the expansion of h . and 77 . 

Results and Discussion of Small-Perturbation Solution 

The formulas derived in the preceding sections of this chapter were programed for 
numerical computation on the IBM 7094 electronic digital computer. The value of y, 
the density ratio across the normal shock, was fixed at a constant value of 0.06 for the 
calculations reported in this and subsequent sections. This choice is justified because 
X varies little with altitude and velocity and the effects of this variation on the stagnation 
solutions are slight. The value x = 0*06 is typical for hypervelocity flight in the atmo- 
sphere of the earth. 

The numerical calculations indicate that the enthalpy is a double-valued function of 
the Dorodnitsyn coordinate 77 in the vicinity of the shock for large values of the Bouguer 
number. An examination of the governing equations failed to show the presence of any 
singularities which might adversely influence the solution in this region when kp is 
large and e small. On the other hand, the results of numerical calculations with varying 
mesh size seemed to rule out the possibility that the doubled-valued behavior can be 
attributed solely to numerical inaccuracies. Consequently, it is suspected that the diffi- 
culty results from truncation of the perturbation expansion and that inclusion of higher 
order terms would either eliminate the problem or increase the value of kp at which 
it first appears. For truncation after the second order term, the condition for validity 
of the solution is ekp < 1. 

Gray gas results .- Shock- layer enthalpy distributions for a gray gas with differing 
values of the radiation cooling parameter e, the Bouguer number kp, the variation with 
enthalpy of the Planck mean mass absorption coefficient Kp (= 9Kp/9h at h = l) and the 
reflectivity of the body surface r w are presented in figures 5 to 7. Although the gray 
gas assumption may not be realistic for most gases of interest, its use is felt to be justi- 
fied in the study of these parameters for two reason s. First, the highly complex and 
v aried spectral structure of absorption -c.Qefficien ts..makes a general p arametric study of 
nongray gases impractical. Second, experience with nongrav calculations indicates that 
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the qualitative dependence of the gray results on the various parameters will c arry o ver 
to mostjiongray cases. 

The decrease in enthalpy level with increasing e is illustrated in figure 5. These 
results indicate that the loss of energy from the shock layer by radiation (that is, radia- 
tion cooling) can produce a noticeable drop in enthalpy for values of e as small as 0.01. 
The dependence of the enthalpy distribution on the Bouguer number (hence, optical thick- 
ness) is also shown in these figures. As expected an increase in the Bouguer number (or 
optical thickness) inhibits shock-layer cooling and leads to higher values of enthalpy near 
the wall. 

The variation of the enthalpy distribution with K p (the enthalpy variation of the 
Planck mean mass absorption coefficient) for several values of the Bouguer number kp 
is shown in figure 6. These effects are most noticeable for optically thin shock layers 
(kp « l.o) and tend to vanish as the optical thickness increases. In a- transparent layer, 
the rate of emission of radiant energy is proportional to the Planck mean mass absorp- 
tion coefficient K p . Thus, gases with small values of Kp (which mean larger values of 
Kp when the nondimensional enthalpy is less than 1) will be cooled more than gases with 
large values of Kp. As the optical thickness increases, smaller K p values still imply 
greater emission rates but they also mean greater absorption and more radiant energy 
available for absorption. The process of absorption tends to counteract the differences 
in emission rates due to differences in Kp. Finally, when radiation equilibrium is 
reached (this state is achieved in the interior of optically thick regions), the energy of the 
particle is independent of its optical properties. Of course, in those regions optically 
close to the shock and the wall, the amount of radiant energy available for absorption is 
not so great as in the interior of the shock layer and particles in these regions cannot 
approach the state of radiation equilibrium (except in a region optically close to a per- 
fectly reflecting surface). Thus, the enthalpy distribution remains dependent on the value 
of Kp near the shock and the wall. This dependence of K p is suppressed near the 
shock where h is almost 1 because the values of K p are nearly the same despite the 
differences in k p . 

The effect of surface reflectivity r w on the shock-layer enthalpy distribution is 
shown in figure 7. If the shock-layer gas is transparent (that is, the gas does not absorb), 
surface reflectivity has no effect on the enthalpy distribution because all photons emitted 
by the layer escape. Whether a photon is absorbed or reflected by the wall is of no con- 
sequence. As the optical thickness of the layer increases, the chance of capture of a 
photon by absorption in the layer is increased. If the surface reflectivity is increased 
also, the probability of capture is increased still further because many photons which 
might otherwise have escaped into the wall are reflected back into the layer and are once 
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again subject to capture there. Consequently, the enthalpy level is higher near a reflecting 
wall than it would be near a nonreflecting wall. 

It can be concluded from these statements that use of a reflecting surface will not 
reduce the radiant heat-transfer rate from the gas to the wall by the factor 1 - r w 
(unless, of course, the gas is transparent) but will reduce it by some smaller fraction. 

This condition exists because the radiant heat flux incident on the wall is larger when the 
wall is reflecting as a result of the higher enthalpy level. In addition, the rate of heat 
transferred to the wall by conduction will be greater, also because of the higher enthalpy 
level. Of course, increasing the surface reflectivity always decreases the total heat- 
transfer rate to the wall because the higher enthalpy level must lead to an increased loss 
of energy by radiation through the shock in the upstream direction and by convection in a 
lateral direction away from the stagnation point. If the energy balance is to be main- 
tained, the rate of heat transferred to the wall must be reduced. 

The effects of variations in the parameters e, kp, Kp, and r w on the rate of 
radiant heat transfer to the wall q^ (normalized by the energy influx to the shock layer, 

1 3\ ' 

are shown in figures 8 to 10. The rate of radiant heat transfer to the stagnation 

LA / 

point was calculated with the formula 

q^ = e(l - r w ) J J** Kp(rj) Bfa) E 2 (k pT^)) d V (127) 

where the optical thickness kpT is given by 

k p T(T7) = k p K p (T}) dr] (128) 


The dashed curves in figure 8 indicate the "no decay limits" for various values of 
the Bouguer number. These limiting curves were computed by assuming the shock 
layer to be isenthalpic so that Kp{rj) = B(ri) = 1. Thus, the no-decay-limit curves are 
given by the formula 


<*5 = K 1 ~ r w) : 


1 - 2E 


s( k p) 


2kt 


(129) 


where E 3 ^kpj is the exponential integral function of third order. This no-decay approxi- 
mation is often used to predict the rate of radiant heat transfer when radiation effects are 

R. 

small. Use of this approximation always gives an upper bound to the true value of q . 

w 

A study of figure 8 indicates that the no-decay limit curve is least accurate in predicting 
the rate of radiant heat transfer in the transparent case kp = 0. This result is expected 
because the enthalpy distribution for the transparent case is the most perturbed from 
an isenthalpic state. Results presented in this figure also indicate the importance of 
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absorption (as characterized by the Bouguer number kp) in reducing the rate of radiant 
heat transfer from the shock layer to the •wall. 

The results presented in figure 9 indicate the differences in k p, the enthalpy varia- 
tion of the Planck mean mass absorption coefficient, are most important when the optical 
thickness of the shock layer is small. 

Here the radiant heat transfer to the wall is greatest for the smallest value of /ip. 
This result, of course, supplements the observation (from fig. 6(a)) that radiation cooling 
is greatest for gases in which Kp is least. The differences in radiant heat transfer to 
the wall brought about by differences in the value of Kp tend to vanish as the optical 
thickness of the layer increases. 

The reduction in radiant heat transfer to the wall due to surface reflectivity is 
shown in figure 10. When the shock layer is transparent, the rate of radiant heat trans- 

O 

f erred q w is in the ratio 1 - r w . However, as the optical thickness of the shock layer 
increases, the ratio becomes somewhat greater than 1 - r w as predicted in an earlier 
discussion of this section. 

The effect of the parameters e, kp, and r w on the shock standoff distance is 
shown in figures 11 and 12. The quantity A is the ratio of the shock standoff distance 
in a radiating shock layer to that in a nonradiating (or adiabatic) shock layer at the same 
flight conditions and was computed with the formula 



The results shown in figures 11 and 12 indicate, as expected, that a decrease in 
enthalpy level (with the consequent increase in density level) in a shock layer leads to a 
reduction in shock standoff distance. 

Nongray results .- The absorption coefficient of high temperature air and other 
gases depends strongly on wavelength. Consequently, the assumption that the gas is gray 
(that is, that the optical properties of the gas are independent of wavelength) is poor indeed, 
and has been resorted to so frequently in the literature only because of the resulting rela- 
tive simplicity. Fortunately, the small perturbation solution derived in this chapter over- 
comes these difficulties by reducing the absorption integrals in the divergence of the 
radiant flux to a form amenable to direct evaluation. Thus, one need only perform an 
integration over a known, albeit complicated, function of wavelength. In view of the cur- 
rent uncertainties, with regard to spectral distributions of gaseous absorption coefficients, 
it was decided to use a simplified model for the absorption coefficient of air. Conse- 
quently, the step function model shown in figure 13 was chosen for use in the calculations. 
The height and width of the steps were chosen so that the simple step function model pro- 
vides an adequate representation of the absorption coefficient of air at a temperature of 
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about 15 000° K as predicted by Nardone et al. (ref. 50) and so that the Planck mean 
absorption coefficients of both distributions are equal. The relative heights of the nine 
steps located at wavelengths less than 0.113 micrometer were chosen to be independent 
of enthalpy whereas the tenth step which covers the wavelength interval (0.113,°°) was 
chosen to vary as the 1.28 power of the enthalpy. The relative heights shown in fig- 
are 13 are for h = 1, where h is the nondime nsional enthalpy. The enthalpy variation 
af the step heights listed is consistent with the condition that the Planck mean mass 
absorption coefficient is proportional to the fourth power of the- enthalpy. 

Shock-layer enthalpy distributions were calculated by using the nongray absorption 
coefficient model for various values of the Bouguer number kp. A comparison of the 
'esults of these calculations with gray calculations by using the Planck mean mass 
ibsorption coefficient is presented in figure 14. 

The maximum monochromatic Bouguer number for the nongray shock layers is 
=.86 times the Planck mean Bouguer number. When the Planck mean Bouguer number kp 
s less than about 0.001 (this case is not shown), the shock layer is optically thin at all 
vavelengths and no perceptible difference between the nongray and the gray calculations 
or the enthalpy distribution can be found. When kp = 0.01, the monochromatic Bouguer 
lumbers for several of the steps are of an order of magnitude of unity and absorption 
tecomes important in the nongray model whereas absorption is still negligible in the gray 
nodel. As a consequence, the enthalpy distribution for the nongray model lies above that 
or the gray model. When kp is increased to 0.1, the disparity between the two solu- 
ions is increased still further. In this case, absorption is very important in those regions 
>f the spectrum for the nongray model in which much of the energy is emitted. Absorp- 
ion is still of minor significance in the gray model. When kp = 1.0, absorption becomes 
mportant in the gray model but still not to the extent that it is in the nongray case. 

Obviously, and not unexpectedly, a gray model which uses the Planck mean mass 
bsorption coefficient will not provide an acceptable estimate of the shock-layer enthalpy 
istribution for a nongray gas unless that gas is optically thin at all wavelengths in which 
significant amount of radiation is transported. Nevertheless, it is very interesting, and 
ncouraging to note that enthalpy distributions computed for the nongray models do not 
iffer significantly in their general shape from those that can be computed for gray models, 
‘hus, it appears that there is some wavelength-averaged absorption coefficient (other than 
le Planck mean when absorption is important but tending toward it in the transparent 
mit) which will provide a good approximation to the enthalpy distribution in a nongray 
as. 

The rate of radiant heat transfer to the stagnation point has been calculated for non- 
ray shock layers. The results are compared in figure 15 with the results of gray cal- 
culations using the Planck mean absorption coefficient. The gray approximation provides 
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a considerable overestimate of the radiant heating even for values of the Planck mean 
Bouguer number as small as 10 " 3 . It is apparent from this result that the tallest steps 
play a very important role in the transfer of energy by radiation. This result is not 
surprising when one considers that nearly 40 percent of the energy emitted by a particle 
in the shock layer is transmitted in the wavelength intervals occupied by the three tallest 
steps. 

It can be concluded from the foregoing discussion that the effective optical thickness 
(or Bouguer number) of a nongray shock layer is greater than that predicted by a gray 
analysis using the Planck mean absorption coefficient. 

OPTICALLY THIN SHOCK LAYERS 
The Transparent Approximation 

Under certain conditions, the Bouguer number, which is indicative of the optical 
depth of the shock layer, is very small compared with unity. When these conditions are 
met, absorption is unimportant and the absorption integrals which are modified by the 
Bouguer number can be dropped from the expression for the divergence of the radiant 
flux vector. (See eq. (87).) This procedure leads to considerable simplification because 
only the local emission rate of radiant energy needs to be considered. All this radiant 
energy is assumed to escape the shock layer and it matters not, insofar as the gas is 
concerned, what path it takes. Consequently, surface reflectivity will have no influence 
on the. enthalpy distribution in the shock layer. Since only the total rate of radiant energy 
emitted locally is of interest, the details of its spectral distribution can be ignored. 

The results of the simplification is the "transparent" form of the divergence of the 
radiant flux vector 

l(V| = -2k p (?7) B(t7) (131) 

The shock layer is termed transparent because the gas is transparent to its own radia- 
tion. Use of the transparent approximation reduces the governing equations from integro- 
differential to a purely differential form. Several investigators (see, for example, refs. 4 
to 7) have taken advantage of this simplicity to obtain approximate analytic solutions. 

The Optically Thin Approximation 

In this paper, a distinction shall be made between the terms "transparent" and 
"optically thin.” A layer of gas will be called transparent if none of the radiation emitted 
by the gas in the layer is reabsorbed. An optically thin layer is one in which a small 
amount of absorption does occur and the optical depth of the layer is small but not zero. 

In the literature, "optically thin" is often used synonymously with "transparent" as defined 
above. 
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P. D. Thomas (ref. 29) expressed concern about the validity of the transparent 
approximation, particularly in the highly cooled region adjacent to the cold wall. The 
transparent approximation is based on the assumption that emission is much greater than 
absorption throughout the shock layer. In regions of small enthalpy, emission no longer 
dominates absorption, and when radiation cooling effects are large, these regions may 
extend over a significant portion of the shock layer. Even when radiation cooling effects 
are small, the value of enthalpy adjacent to the wall tends to vanish^ and absorption must 
become important compared with the local rate of emission. Of course, for this case, 
the region of nonvalidity is very small and has no appreciable effect on overall proper- 
ties such as the radiant energy flux to the wall and the shock standoff distance. 

Thomas sought to modify the transparent equations in order to take into account 
this small amount of reabsorption. He did so by expanding the Planck function B^(t) 
which appears in the integrand of the divergence of the radiant flux vector in a Taylor 
series about the zero of the argument of the displacement kernel Ej^kpj - t^|^. The 
expansion is then arbitrarily truncated after the linear term. • Strictly speaking, this 
procedure can be used only when the Planck function varies slowly within a photon mean 
free path length. Obviously, this criterion is not met when the shock layer is optically 
thin (particularly close to the wall, the region of greatest interest, where the enthalpy and 
hence, the Planck function varies rapidly) and some doubt must be cast on the validity of 
Thomas' analysis. 

It would appear that the effects of small absorption could better be discovered 
through a straightforward expansion of the equations in terms of the Bouguer number kp. 
Such a solution, up to first order in kp, is presented here. In order to simplify the 
analysis, the exponential integral functions E£(x) and E 3 (x) which appear in the expres- 
sion for the radiant flux are replaced by the exponential functions e"^x an( j (l/2)e~^ x , 
respectively. The particular form of the exponential functions was chosen so that the area 
under the curve of E£(x) and the approximating exponential are equal for the interval 
(o, x a)> f° r x a « 1; thus the value for the radiant flux reduces to the proper value in the 
transparent limit. This substitute kernel approximation has been used with considerable 
success in a variety of problems of radiant transfer. (See, for example, refs. 11, 32, 33, 
and 51.) 

Use of the substitute kernel approximation reduces the expression for the radiant 
. flux to the form 


I ^An element of gas approaching the wall requires an infinite time to reach its 
.(destination. Because of this time required and the fact that the rate of energy lost by 
d radiation is proportional to a positive power of the enthalpy, the enthalpy of a transparent 
a gas must approach zero as the particle approaches the wall. 
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The divergence of the radiant flux vector is 
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Here the monochromatic optical depth is 

rV 

r x (v) = J K X \V) d?7 


( 133 ) 


( 134 ) 


It was seen in the section "The Small-Perturbation Solution" that an expansion of 
the governing equations in terms of the small parameter e led to a fortuitous uncoupling 
of the energy and momentum equations. Unfortunately, the same is not accomplished when 
the expansion is performed in terms of kp. It is frequently pointed out in the literature 
(for example, refs. 4 and 6) that the coupling is very weak. The solution to the momentum 
equation depends on the solution to the energy equation through the term a ^3(77) (see 
eq. (72)) where a 2 is order of magnitude x* An analysis of equation (72) indicates that 
the contribution of this term to f (77) is order \fx- Since x I s very small (typically 
0 . 06 ), the effect of the solution of the energy equation to f (77) is order 0.25. Advantage 
can be taken of this situation by replacing ^3(7?) = h (??) which appears in the momentum 
equation (72) with h, the integrated average of h (77) over the interval (0,77^; that is, 

H = h( 7 ?) d ?7 ( 135 ) 

^A J 0 

(This replacement is tantamount to solving the momentum equation by the integral method 
with f (77) assumed to be a quadratic function.) This approximation has the twin virtues 
of retaining the coupling, albeit in approximate form, and greatly simplifying the solution 
to the energy equation. 
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Now the governing system of equations takes the approximate form: 


f(? 7 ) h'(r?) + el[? 7 ] = 0 

(136) 

2f(?7) ~ [f't 7 ?)] + = 0 

(137) 

o 

1! 

& 

(138) 



(139) 
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(140) 


h(7? A ) = 1 (141) 

where l£Vj is given by equation (133). When the Bouguer number is very small, the 
absorption integrals in equation (133) assume a secondary significance throughout the 
domain of the problem. Neglecting these absorption integrals reduces the system to a 
purely differential form. If, as expected, when the Bouguer number kp is small, the 
presence of absorption only slightly influences the solution, one can, to reasonable accu- 
racy, evaluate the absorption integrals by using the transparent solution for h. The 
perturbation expansion scheme used herein follows the general outline discussed above. 
Mathematical details are presented in appendix B. 

The zero-order, or transparent, solution is 


where 
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(147) 


fio-1 V x > dx 


'0 

It has been shown that the Planck mean mass absorption coefficient normalized by 
its value immediately behind the shock can be adequately represented by 


K-p(h) = 


fh^ 

IIV 

[ch yi 

(h < h*)_ 


(148) 


1 2 

where - W^h* is the value of the enthalpy (depending on the pressure, of course) at 
which the value of the exponent of h changes. The constant C is obtained by equating 
the two expressions for «: p (h) at h = h* with the result 


C = (h*) 


72 ~n 


(149) 


It also has been shown that the nondimensional Planck function B(h) is approximately 
given by the expression 

B(h) = h 2 ' 2 <150) 


When the correlation formulas (148) and (150) are introduced into equation (142), 
the integration on the left-hand side can be carried out, and the solution for hg(x) given 
by the explicit formula 
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h Q (x) 
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for h 0 (x)gh*,and 


h Q ( x ) 
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for hQ(x) < h* . 

The first-order solutions which include the effects of absorption, surface reflec- 
tivity, and nongray radiation are: 
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The Poincare-Lighthill-Kuo Solution 

Careful inspection of the last term on the right-hand side of equation (152) reveals 
that the first-order term h^(x) displays a singular behavior near rj = 0 (where ho 
approaches zero). By way of illustration, consider the case of a gray gas with 
Kp(h) = h y and B(h) = h 2 - 2 . In this case the term in question is proportional to the 
quantity 


f 1 h^ +2 * 2 r 

v+2.2 _d|_ = ^0 r h -(r + 3.4) . J 

0 Jh 0 hr+ 4 - 4 y + 3.4L° J 


Near the wall, h Q approaches zero and equation (156) approaches — - which 

y + 

increases without limit. This seemingly anomalous behavior can be explained as follows: 
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The first-order solution represents a gas which absorbs radiation at a rate determined 
by the absorption coefficient for the zero-order solution (the magnitude of the incident 
radiation is independent of the amount absorbed) while it emits energy at a rate propor- 
tional to the derivative with respect to hg of the emission rate for the zero-order solu- 
tion. Both the absorption and emission rates tend to zero as an element of gas approaches 
the wall. However, the emission rate tends to zero much more rapidly than the absorption 
rate. The difference in the limiting behavior of these rates coupled with the infinite resi- 
dence time for an element of gas in the stagnation region allows the gas element to absorb 
an infinite amount of energy and thus the enthalpy of the gas adjacent to the wall becomes 
infinite . 

The difficulty which has arisen as the result of the singularity can be avoided 
through the use of the P-L-K perturbation of coordinates procedure which transforms the 
coordinate in such a way that the singularity is removed from the boundary (at x = 0) to 
a point outside the domain of the problem (a slightly negative value of x) . Mathematical 
details of the application of this method are described in appendix B. The P-L-K solu- 
tions to first order in kp are 

x = y + k p xj(y) (157) 

h ( x;k p) = h 0< y > U58) 

f(x; k p) = fj(y) + k p fj(y) (159) 

where the starred coefficients in the P-L-K expansions are related to the unstarred coef- 
ficients in the regular perturbation expansions. (See ref. 49.) 

hj(y) = h Q (y) (160) 

f£(y) = f 0 (y) (161) 

fj(y) = fj(y) + Xj(y) fg(y) (162) 



-h 1 (y) 

h o (y) 


(163) 


Results and Discussion of Solution for Optically Thin Shock Layers 

In order to obtain some indication of the accuracy of the optically thin shock-layer 
approximation, the results computed for a typical case are compared in figure 16 with 
the results computed by means of the small-perturbation method and the results of a 
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numerical calculation performed by Howe and Viegas (ref. 9). The agreement 
the three solutions is excellent. However, a word of caution should be interject 
in order to avoid the implication that the numerical results of Howe and Viegas r 
sent the "exact" solution to the inviscid, plane-parallel geometry, stagnation flow . 

The results of Howe and Viegas include viscosity, heat conductivity, and body curva 
ture. The effects of curvature are expected to be very small. The flight conditions 
(Woo = 9.75 km/sec; p s = 10 atm) were chosen to insure that the boundary layer was very 
thin so that "displacement" effects on the inviscid region were minimized. Finally, the 
thermodynamic and optical properties used by Howe and Viegas were obtained from their 
own correlations while the optically thin and small perturbation methods were computed 
by using the correlations presented herein. Thus, the comparisons between the results 
from the methods of this paper and the numerical results of Howe and Viegas are as 
much, or more, checks on the validity of using the inviscid approach and checks on the 
similarity of two different sets of correlations as they are checks on the accuracy of the 
analytical methods of this paper. It is conceivable that errors due to the various factors 
mentioned tend to cancel in this example. Nevertheless, the individual errors due to the 
omission of viscosity, heat conductivity, and curvature and due to the difference in cor- 
relation functions are expected to be very small so that the excellent agreement can still 
be interpreted as an indication of the accuracy of the small perturbation and optically thin 
methods . 


The optically thin solution was used to study the effects of the radiation cooling 
parameter e, the Bouguer number kp, the surface reflectivity r w , and the enthalpy 
dependence of the absorption coefficient on the shock-layer enthalpy distribution, the rate 
of radiant heat transfer to the stagnation point, and the shock standoff distance. As in the 
previous section, the density ratio y across the near-normal section of the shock was 
fixed at a value of 0.06. In addition, all the results are limited to the case of a gray gas. 

The effect of absorption on the enthalpy distribution is indicated by the curves of 
figure 17. The solid curves represent the enthalpy distributions in a transparent shock 
layer for e = 0.01, 1.0, and 100. The dashed curves represent the enthalpy distributions 
: in optically thin shock layers for the same values of the radiation cooling parameter e. 
Values of the optical thicknesses are shown in the figure. These results show the expected 
trend, the enthalpy level falling as the radiation cooling parameter e increases. Absorp- 
, tion tends to increase the enthalpy particularly in the cooler regions of the flow. Absorp- 
! tion also affects the location of the shock and reduces the value of 77 . (the location of 
the shock in terms of the Dorodnitsyn coordinate) because of the decreased density level. 
Although the value of 77 ^ decreases, the shock standoff distance A increases with 
I increasing optical depth. 
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The effect of the enthalpy dependence of the absorption coefficient on the enthalpy 
distribution in transparent shock layers (kp = o) is shown in figure 18. In figure 18(a) 
the absorption coefficient was given by the relation Kp = h^, where y takes on the 
values^ 3, 4, and 5. The value of y determines how the rate of energy emission varies 
with enthalpy across the shock layer. The rate of energy loss by radiation will decrease 
more rapidly as the enthalpy falls if y is large than if it is small. Consequently, the 
enthalpy distribution for a large value of y lies above that for a smaller value. This 
trend, of course, is the same trend exhibited by the small-perturbation solutions of the 
preceding section. In figure 18(b) the absorption coefficient is given by the relation 

Kp = Ch y where C = (h ) and y = y^ = 4 for h < h* , and C = 1 and y = y^ = -1 

for h > h* . This model should be used when the shock-layer temperatures are in excess 
of about 20 000° K since at moderate altitudes n (the enthalpy at which the exponent y 
changes value) corresponds to temperatures of approximately this value. The effects of 
varying h* are shown in figure 18(b). A decrease in h* produces a decrease in the 
enthalpy level because y takes on the smaller value (-1) throughout a greater portion of 
the shock layer. 

The effect of surface reflectivity r w on the enthalpy distribution is shown in fig- 
ure 19. Of course, this effect vanishes in a truly transparent layer. With a small amount 
of absorption an increase in reflectivity brings about an increase in enthalpy level, the 
greatest increases occurring adjacent to the wall. These results corroborate the findings 
obtained with the small -perturbation solution. 

The variation with the radiation cooling parameter e of the rate of radiant heat 
transfer to the wall for various values of the Bouguer number is presented in figure 20. 
Also shown in this figure are two limit curves. One of these curves is labeled the "no 
decay limit" and was computed by assuming that the shock layer was isenthalpic and 
transparent. The second limit curve is labeled the "available energy limit” because it 
represents an upper bound to the radiant flux on the basis of energy balance. The amount 
of energy entering the shock layer per unit time per unit area of the shock surface has 
been normalized to unity. If all this energy is radiated out of a transparent shock layer, 
only one -half will be incident on the wall. 

The curve labeled kp = 0 shows the effect of radiation decay or cooling in reducing 
the rate of radiant heat transfer to the wall. The remaining curves indicate the important 
effect of absorption (as characterized here by the Bouguer number kp) in reducing the 
rate of radiant heat transfer to the wall. Although values of kp presented in figure 20 
are as large as 3, the corresponding shock layers are all optically thin « lj. 


®These values of y are typical for air at temperatures less than about 20 000° K._ 
(See, for example, eq. (94).) 


The effect of the enthalpy dependence of the absorption coefficient on the rate of 
radiant heat transfer to the wall in a transparent shock layer (kp = o) is shown in fig- 
ure 21. It is apparent that an increase in the exponent y (which appears in the correla- 
tion formula Kp = h^) magnifies the effect of decay on the rate of radiant heat transfer. 

The effect of the radiation cooling parameter e on the shock standoff distance for 
various values of the Bouguer number kp and y is shown in figure 22. As expected, 
an increase in e reduces the value of A (the ratio of shock standoff distances with and 
without radiation) for given kp and y because the cooling by radiation tends to increase 
the density level in the shock layer. Increases in kp and y for fixed e inhibits the 
effect of decay on A whereas these increases magnified the effects of decay on the rate 
of radiant heat transfer. 

The variation of shock-layer optical thickness kpT^ with the radiation cooling 
parameter e and the Bouguer number kp is shown in figure 23. When the absorption 
coefficient varies as a positive power of the enthalpy, the shock-layer optical thickness 
may be very much less than one even if the Bouguer number is order of magnitude one or 
greater provided that e is sufficiently large. 

The criterion for the validity of the analysis presented in this chapter is that the 
optical depth of the shock layer be much smaller than one for those wavelength regions 
in which a significant amount of radiant energy is transported. It has been shown herein, 
for the case of a gray gas for which the absorption coefficient is proportional to a posi- 
tive power of the enthalpy, that this condition is always less restrictive than the condition 
that the Bouguer number kp is very much less than one. However, for the more realis- 
tic case of a nongray gas, the criterion stated is generally more restrictive than the con- 
dition kp « 1. In mathematical terms the criterion implies the inequality 

2(1 + r w )k p k x (77) k x U) B X U) d^j dX « 1 (164) 

The quantity on the left-hand side of the inequality is the first-order term in the expansion 
of ijjfj, the divergence of the radiant flux vector, in terms of the Bouguer number kp. 
When both k x and B^ are proportional to a positive power of the enthalpy, an upper 
bound to the aforementioned quantity can be obtained by replacing K x (rj), k x (£), and 
B x (£) by their values at V = immediately behind the shock. The result is 

2(1 + r w )k p 77 A 4( V A ) B x (t? a ) dX « 1 (165) 

If the same substitution is used for a gray gas, the result is simply 

(166) 


2(! + r w )k p ?7 A « 1 
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because k ^\ 

step function model for the absorption coefficient of air, which was used in the small 
perturbation solution (see fig. 13), is used to evaluate the quantity (165), the result is 
about 60 times greater than the corresponding gray quantity given by equation (166). 

Thus, the criterion for the validity of the optically thin analysis, in this nongray example, 
is 


/-» OO 

^ and j ^ are both identically equal to 1. When the nongray 


60k p « 1 


for small values of the radiation cooling parameter e. For larger values the criterion 
could probably be relaxed somewhat ^for example, 60kpT^ « 1^. As a result of this 
criterion, the practical applicability of the optically thin analysis (and consequently of all 
transparent analyses) is seriously restricted. 


THE OPTICALLY THICK SHOCK LAYER 


The Optically Thick Approximation 

A qualitative description of the optically thick shock layer has been given by 
Goulard (ref. 5). He pointed out that this layer is characterized by an isothermal region 
between two thin boundary layers adjacent to the shock and the wall. The boundary layer 
immediately behind the shock is a result of the cooling of the hot gas by radiation through 
the transparent shock. Because radiation travels only a short distance before being 
absorbed in an optically thick layer, this energy loss is restricted to a narrow region 
which extends approximately a photon mean free path. Once this initial adjustment in 
energy has occurred, the gas particle is carried into the interior of the shock layer by the 
flow where convection is the dominant mode of energy transport. In this region, the 
enthalpy of the gas is essentially constant. As the particle nears the cold wall, moving 
ever more slowly as it does so, convection becomes of decreasing importance and energy 
transfer by radiation begins to assume the major role. Finally, in the immediate vicinity 
of the wall, all the energy transport proceeds by means of radiation. When the emissive 
power in the interior (or isothermal portion) of the shock layer is large, the "take-over" 
by radiation occurs at greater distances from the wall than if the emissive power is small. 
Thus, the thickness of the wall boundary layer depends not only on the optical thickness of 
the shock layer but also on the emissive power of the gas. 

Although the shock layer is optically thick, the thermal boundary layer immediately 
behind the shock is not and thus the Rosseland or diffusion approximation so commonly 
used in the study of optically thick gases cannot be applied in this region. The Rosseland 
approximation is valid only in regions of an opaque gas which are at great optical distances 
from all radiation boundaries (a perfectly reflecting barrier is not a radiation boundary) 
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and in which the thermodynamic and optical properties do not vary greatly within a photon 
mean free path. Neither of these conditions are met in the shock boundary layer. 

The conditions of validity for the Rosseland approximation might hold throughout 
much of the wall boundary layer if the emissive power of the gas is sufficiently large. 
However, the approximation must break down optically close to the wall. The use of a 
temperature jump boundary condition as suggested by several investigators (refs. 27, 52, 
and 53) has proven successful in problems of radiant and combined radiant and conductive 
energy transport. Whether this concept can be applied with equal success to problems of 
combined radiant and convective energy transport has not, as yet, been demonstrated. In 
a region optically close to a radiation boundary, the temperature predicted through the 
use of the Rosseland approximation and a slip boundary condition represents not the tem- 
perature of the molecules of the gas, but a sort of average photon temperature. The con- 
vective heat flux depends on the molecular temperature. Thus, it is not clear that the 
slip boundary condition can be used in a problem of combined radiant and convective 
energy transport. There is a basis for optimism when considering the problem of this 
chapter, however, in that the convective flux may be negligible compared with the radiant 
flux optically close to the wall. 

In order to arrive at a solution to the problem of the optically thick shock layer, the 
substitute kernel approximation, introduced earlier is used. It is shown that in the 
interior of the shock layer and close to the wall but not in the shock boundary layer, this 
method is equivalent to using the Rosseland approximation with slip-boundary conditions. 
The use of this approximation will restrict the analysis to gray gases. This restriction 
is not a condition for application of these approximations, but has been invoked here to 
avoid the considerable additional complication that relaxation of this restriction would 
incur. 

The Substitute Kernel Approximation 

In this and the subsequent section, it will be convenient to rewrite the energy equa- 
tion (eq. (71)) with the optical path length r as the independent variable; that is, 

f ( t) h'(T) + eljV] = 0 (167) 

Here, f is the nondimens ional stream function; h, the nondimensional enthalpy; and e, 
the radiation cooling parameter. The divergence of the radiant flux vector l|~rj is 
given by the expression 

l[r] = l[rfj = -2 B(t) E 2 (0) + k p Jj* B(t) E^kpjr - t|) dt + 2k p r w E 2 (k p r) B(t) E 2 (k p t) dt (168) 
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T is the nondimensional absorption coefficient; B, the nondimensional Planck 
ody function; kp, the Bouguer number; r w , the reflectivity of the wall (at r = 0); 
value of the optical path length at the shock; and Ej and E 2 , the exponential 
llll _i functions of first and second order, respectively. 

In order to simplify the analysis, the substitute kernel approximation is used. For 
the optically thick shock layer, the appropriate substitution for E 2 (x) is found to be 
(3/4)e“^ // ^ . This substitution satisfies the conditions that the areas under the two func- 
tions over the domain 0 = x S °° are equal and that the expressions for the radiant flux 
approach the Rosseland expression as x increases without limit. 

If the expression for the radiation flux is differentiated twice with respect to r, 
the integral terms can be eliminated with the result 

i ”[ T ]"! k P i [ r ] = -I B " (r) (i69) 

The energy equation (167) can then be used to eliminate ijVJ 

[f(r) h’(r)] ” - | k|f(r) h ’( t) = | eB”(r) (170) 


This alternate form of the energy equation is a third-order nonlinear ordinary differential 
equation, the solution of which must satisfy the condition h(l) = 1. Two additional con- 
stants of integration are introduced by the solution of equation (170). These constants 
are determined by satisfying appropriate physical conditions or by satisfying the original 
integro -differential equation (eq. (167)). 


An expression for the flux of radiant energy which enters the wall can be obtained 
very simply. The expression for the flux incident on the wall is 




• — k t 
2 P 


dt 


When the integro-differential form of the energy equation (eq. (167)) is evaluated at 
it becomes (since f(0) =0 is a boundary condition) 


B(0) - | k p (l + r w ) J Ta B(t) e ^ dt = 0 


r = 0, 


Thus, the flux entering the wall can be written in terms of the value of the black-body 
emissive power of the gas adjacent to the wall; that is, 



1 - r. 


w 


1 + r 


w 



(171) 
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Boundary-Layer Analysis 

In terms of the substitute kernel approximation, the complete differential system 
governing the flow in the stagnation region of a radiating shock layer is 


h’(ry] " - | eB”(r) - | k|f(r) h'(r) = 0 

(172) 

2f (V) f"(7) - |/(7)] 2 + a 2 h(7?) = 0 

(173) 

f(0) = 0 

(174) 

f W Bl 

(175) 

- - a . 
' ' fox(l - x) 

(176) 

»{ t a) = 1 

(177) 


f(r) h'(-r) 






k P r w e 
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(178) 


When the optical thickness of the shock layer is such that 
equation (172) becomes asymptotic to the equation 

f(r) h'(r) =0 


kp » 1 


and kl » e, 


(179) 


Thus, the enthalpy approaches a constant. It can be shown by attempting to satisfy equa- 
tion (178) as a condition, that this solution is valid only at large optical distances from 
both the shock and the wall (unless it is a perfectly reflecting wall). It also becomes 
clear that the value of this constant, hereafter denoted h2, cannot be determined without 
knowledge of the shock boundary layer. 


If the enthalpy throughout most of the shock layer is constant with a value h 2 , the 
density will be constant also with a value pg. In this case, the momentum equation may 
be easily solved with the result that 


f(p) = 


- ap A 



+ ap A 


,l£2 /J2_\ 

Hvw 


(180) 


53 



The first approximation to the shock standoff distance is 

( p s/ p 2> h 2 X - (181) 

Rs 1 + \j 2 (Ps/P 2)xp " X) 1 +^/2h 2 x(l - x) 

In addition to the region of constant enthalpy in the interior of the layer, there are 
two thermal boundary layers, one immediately behind the shock and the other adjacent to 
the wall. 

The forms that the energy equation assumes in these boundary layers can be deter- 
mined by means of conventional boundary-layer techniques. In the vicinity of the shock 
the "stretched" coordinate 

t = (r A - r)k p (182) 

is introduced. Close to the shock the quantity f(T) is slowly varying and may be ade- 
quately represented by the first term in the Taylor expansion about £ = 0; that is, 


f(£) ~ 1 (^a) = i (183) 

Substitution of equation (182) and equation (183) into the energy equation (172) gives the 
shock boundary-layer equation 


h"(S) + 1 B ’ ( ^ " I hU) = Constant (184) 

Solution of this equation is complicated by the nonlinear term B'(£). If e is at 

Kp 

least an order of magnitude less than kp, this term can be neglected and the solution 
to equation (184) is easily found. This solution is 


m = (i 


- h 2 y 


■k 

2 ? 


+ h 2 


(185) 


The constant h 2 can be determined by writing an energy balance across the shock 
boundary layer. This energy balance is 


1 =h 2 


3 _e_ r 

4 kp Jo 


- 3 | 

BU) e 2 d^ 


(186) 


When e « kp, condition (186) reduces to 


h 


2 


* 1 


2 kp 


and it is apparent that h 2 approaches one and the boundary layer ceases to exist. Thus, 
there cannot be a shock boundary layer with a thickness characterized solely by the optical 
path length in the gas. 
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An approximate solution to the boundary-layer equation (eq. (184)) can be obtained 


if the nonlinear term 




is replaced by an appropriate linear term; for example, 


! £ B ’ (4> ’ f £ + Bh(4 >] - 1 £ Bh '< 4 > 

where the constant B is arbitrary and represents a mean value of the black-body emis- 
sive power B and the constant B is arbitrary and represents a mean variation of B 
with h over the range of values of h encountered in the shock boundary layer. The 
linearized version of equation (184) has the simple solution 


h(£) = (l - h 2 )e 1 


+ h 2 


where 




l + l + 


b 


Where e/kp is very much less than unity, this solution reduces to equation (185), 
e/kp is very much larger than unity, the solution takes the form 

3e -B| 


h(£) = (l - h 2 )e 2kp + h 2 


(187) 


(188) 


When 


(189) 


and, in this limit, the thickness of the shock boundary layer is determined by the param- 
eter e ~ * instead of simply kp . Thus, the shock boundary layer can be very much 
thinner than a photon mean free path if the black-body radiative power behind the shock is 
large. This effect was shown by Heaslet and Baldwin (ref. 33) in their study of radiation- 
resisted shock waves. Simply stated it means that a particle starting immediately behind 
the shock loses energy at such a rapid rate by means of radiation that it is substantially 
cooled in the time that it takes to travel only a small section of a photon mean free path. 

A value for the constant B can be obtained from the condition 

"i-|£( l + h *)“i-|-° < 190) 

This condition was derived by integration of the nonlinear energy equation (eq. (184)) 
between the limits zero and infinity and substitution into the result of the linearized solu- 
tion (eq. (187)). In addition, the correlation formula B = h 5 was used. It was shown 
earlier that 6 ~ 2.2. However, the ensuing analysis will be greatly simplified, without 
any significant loss in accuracy, by setting 6 = 2. The value of B so obtained is 


B= 1 +h 2 


(191) 
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A second condition is required to determine uniquely the enthalpy distribution in 
the shock layer. The energy balance relation (eq. (186)) evaluated with the aid of the 
linearized solution provides this condition, which is 


S >2 - 2 ^ - h 2 ) + |(! + 21 > 2 ) - 4 ^(1 - h 2 ) 


<JL>1 + 




= o 


(192) 


The quantity can be eliminated between the conditions (190) and (191) to obtain an 
expression for h 2 , the enthalpy level in the interior of the shock layer, as a function 
of e /kp. The result of this calculation is presented in figure 24. 

The thickness of the shock boundary layer (in terms of optical path length) is char- 
acterized by the parameter ^w^kp^ - '*'. A plot of coj as a function of e /kp is pre- 
sented in figure 25. 

As has been indicated previously, there is also a thermal boundary layer due to 
radiation adjacent to the wall. If this boundary layer is thin, which shall be assumed, 
herein, the dimensionless stream function f(r) may be represented by the first two 
terms of its MacLaurin expansion 

f(r) * f(0) + rf’(0) 


By employing the kinematic boundary condition (174) and the asymptotic solution for f ( 77 ) 
(eq. (180)), one finds 


f(r) 


k( 0 )\/h2 


r = br 


(193) 


Of the several approximations introduced in the analysis of this section, this approxima- 
tion is perhaps the poorest because the requirement that the wall boundary layer be thin 
with respect to the optical path length t does not necessarily imply that it is thin with 
respect to either the Dorodnitsyn coordinate 77 or the geometric coordinate z. 


Substituting this expression into the energy equation (eq. (172)), introducing the 
"stretched" coordinate 

? = Tk* /2 (194) 

and neglecting terms of order kp* yields the boundary-layer equation 


k T 


B”(0 +|-f b?h'(?) = 0 


(195) 


In general, equation (195) is nonlinear and does not possess an analytic solution. A sim- 
ple approximate analytic solution can be obtained by replacing the quantity h’(£) with 
hB’(£), where h is an as yet undetermined constant. This substitution reduces equa- 
tion (195) to the linearized form 
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where 


B"(?) + 2w 2 ?B’(0 = 0 


(196) 


<2 4 e 


(197) 


The solution to equation (196) is easily found with the result 

B(0 = B w + (B 2 - B w )erf (fJ^ (198) 

The quantities B w , the nondimensional black-body emissive power of the gas 
adjacent to the wall, and (because it contains the arbitrary constant h) are still 

unknown. One condition for evaluating these quantities can be obtained by integrating the 
nonlinear boundary- layer equation (195) with respect to £ between the limits zero and 
infinity. In performing this integration, it is convenient to eliminate the term B"(£) in 
equation (195) with equation (196). Then it is found that 


g h 2 - h w _ i 
B 2 “ B w h 2 + h w 


(199) 


Here B 2 is the nondimensional black-body emissive power of the gas in the interior of 
the shock layer and h w is the nondimensional enthalpy of the gas adjacent to the wall. 
The last equality in expression (199) holds because it has been assumed that B = h^. 


The second condition is obtained by evaluating the integral condition (eq. (178)) at 
r = 0. When written in terms of the boundary-layer coordinate this condition becomes 


B„ = 2 k l/2 


4 P 


p CO 

( 1 + r w ) J 0 B (0 


-H /2? 

e 2 df 


( 200 ) 


Substitution of the linearized solution for B(C) into equation (200) and integration yields 


(! + r w ) 


9kp 

16 ct)c 


B w = 


erfc\ 


9kj 
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16w, 
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(l - r w ) + (l + r w )e 1 erfc 
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( 201 ) 


Equations (197) and (199) can be used to eliminate oj 2 from equation (201) to yield a 
transcendental equation for the value of the black-body emissive power at the wall B w . 

The variation of B w as a function of the radiation cooling parameter to Bouguer 
number ratio e/kp for various values of kp and the exponent y (from the correlation 
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formula k = h^) is shown in figure 26. This curve has particular significance because 
the ratio of radiant heat transfer to the cold wall is directly related to B w through 
expression (171). The variation of the quantity cog (eq. (197)) with these same param- 
eters is presented in figure 27. 

The Rosseland Approximation 

The Rosseland or diffusion approximation is frequently used in the study of prob- 
lems in which the medium is optically thick. As was pointed out earlier, this approxi- 
mation is not valid in regions optically close to a radiation boundary or in regions in 
which the optical and thermodynamic properties vary significantly within an optical path 
length. Some investigators have attempted to correct the former deficiency through the 
use of temperature-jump boundary conditions and have achieved considerable success in 
problems of pure radiant or combined radiant and conductive energy transport. 

In this section, an attempt is made to use the Rosseland approximation and 
temperature-jump boundary conditions to analyze the optically thick shock layer. It is 
hoped that this exercise will provide some insight into the attributes and limitations of 
this approximation in problems of combined radiant and convective energy transport. 

With the Rosseland approximation for the divergence of the radiant flux, the energy 
equation becomes 

k^ 

B"(t) + | ~ f(r) h’(r) = 0 (202) 

This equation is the same as equation (172) except for the omission of the third-order 
differential term jf(r) h’(r)J . 

In the interior of the optically thick shock layer, equation (202) reduces to 

f(r) h'(r) = 0 

provided kp » e/kp. This result is identical to the result obtained by means of the 
substitute kernel approximation. This agreement is not surprising because the diffusion 
approximation is known to be valid in this region. Of course, the value of the constant 
enthalpy in the interior of the shock layer cannot be determined until something is known 
about the shock boundary layer. 

If the usual type of boundary-layer analysis is applied to the energy equation in the 
Rosseland approximation (202) for the neighborhood of the wall, the result is 

B"(0 +|^b?h'(?) = 0 
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This equation is identical to the wall boundary-layer equation in the substitute kernel 
analysis. Two boundary conditions are required to determine the solution to this equa- 
tion completely. One of these conditions is 


lim B(0 = B 2 

Here B 2 is the black-body emissive power of the gas in the interior of the shock layer. 
The second is the jump boundary condition, written in terms of the black-body emissive 
power B rather than the temperature 

B w = CB’(O) = §(^Cq* (203) 

I 

The second equality follows from the expression for the radiant flux in the Rosseland 
approximation. The constant C is usually evaluated by requiring the flux to be correct 
in the black-body limit. (See, for example, ref. 37.) However, it is noted that condi- 
tion (203) is identical to the condition used in the substitute kernel approximation (that 
is, eq. (200)) if C is chosen to be 


2 1 + r w 

3kp 1 - r w 


(204) 


Thus, the results obtained in the wall boundary layer by the two methods are identical if 
C is chosen to satisfy equation (204) . It has been shown that the two methods also give 
identical results when applied to the problem of combined radiative and conductive energy 
transport between two infinite parallel plates separated by a radiating and conducting gas. 
(See ref. 37.) 


If the usual boundary-layer analysis is used to obtain the boundary-layer form of 
the energy equation in the Rosseland approximation for the neighborhood of the shock, the 
result is 



= Constant 


(205) 


This equation is not identical with the shock boundary-layer equation in the substitute ker- 
nel approximation because of the omission of the third-order differential term. Inspec- 
tion of equation (205) indicates that any solution other than the trivial solution h(|) = h 2 
will not tend to a constant h 2 as £ becomes very large. Thus equation (205) cannot 
be forced to satisfy simultaneously the conditions h(0) = 1 and lim h(|) = h 2 . Appar- 

ently then, the jump boundary condition at the shock must be h(0) = h 2 , but this result 
leads nowhere as there is insufficient information to determine h 2 accurately. 

It must be concluded then that the Rosseland approximation with slip -boundary con- 
ditions is not sufficient by itself to be used in the analysis of the complete shock layer. 
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It can be used in the combined radiation and conduction problem because the two separate 
energy fluxes are represented by similar mathematical models and may be treated as an 
equivalent radiation- alone or conduction- alone problem. Even in this case, it is not pos- 
sible to obtain a temperature distribution nor to determine separately the radiant and 
conductive fluxes optically close to a boundary. This inability to determine a tempera- 
ture distribution optically close to a boundary (such as a transparent shock) presents a 
serious obstacle to the solution of the combined radiation and convection problem because 
the convective flux depends on the temperature distribution. 

Results and Discussion of Solution for Optically Thick Shock Layers 

The analysis presented in this section applies only when the two enthalpy boundary 
layers are thin compared with the shock standoff distance with these distances expressed 
in terms of the Dorodnitsyn variable rj. When e/kp, the ratio of the radiation cooling 
parameter and the Bouguer number, is much less than one, the thickness of the shock 
boundary layer is characterized by the inverse of the Bouguer number k^; whereas 
when e/kp is large, the shock boundary-layer thickness is characterized by the inverse 
of the radiation cooling parameter e~*. The thickness of the wall boundary layer is char- 

/ / 2\l/2 

acterized by the parameter I e/kp I . Thus, the most restrictive conditions on the 

applicability of the optically thick analysis are that kp » 1 for e small and kp » e*/ 2 
for e large. 

Several enthalpy distributions were calculated with the formulas developed in the 
preceding section. The results are presented in figure 28. The previous discussion of 
the effects of the parameters on the thicknesses of the boundary layers is substantiated 
by these results. The effect of the Bouguer number kp and the ratio of radiation 
cooling parameter to Bouguer number e/kp on the shock-layer optical thickness kpr^ 
is also shown. The effect of e/kp depends on the enthalpy variation of the absorption 
coefficient. In the cases shown, the absorption coefficient is proportional to the fourth 
power , of the enthalpy and an increase in e/kp brings about a reduction in the shock- 
layer optical depth. The nondimensional value of the enthalpy of the gas adjacent to the 
wall (which is related to the rate of radiant heat transfer to the wall through eq. (171) 
and the correlation formula B = h^) decreases with increasing e/kp and/or kp. 

The effect of y (where y is the exponent in the correlation formula k = h’ 5 ') and 
the surface reflectivity r w on the character of the wall boundary layer has not been 
shown but can be readily deduced from the curves of figures 26 and 27. Increasing y 
tends to reduce the optical thickness of the wall boundary layer and increase the value of 
the enthalpy of the gas adjacent to the wall. It can be shown that the wall boundary-layer 
thickness expressed in terms of the Dorodnitsyn coordinate T} is only slightly affected 
by a change in y. Increasing the surface reflectivity r w tends to increase the optical 
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thickness of the wall boundary layer and increase the value of the enthalpy of the gas 
adjacent to the wall. When expressed in terms of the Dorodnitsyn coordinate rj, the 
boundary-layer thickness decreases with increasing r w . These results are consistent 
with the results obtained with the small-perturbation solution. 

The manner in which the rate of radiant heat transfer to the wall q depends on 

w 

the radiation cooling parameter to Bouguer number ratio e/kp, the Bouguer number kp, 
the variation with enthalpy of the absorption coefficient (through the exponent y of the 
correlation formula k = hT), and the surface reflectivity r w is indicated in figure 26. 

"R 

For fixed values of kp, y and r w , the rate of radiant heat transfer to the wall q 

increases with increasing e /kp. It appears that q w would become asymptotic to the 

available energy limit of 1/2 as e/kp increased without limit. As the Bouguer number 

kp increases (hence, the shock-layer optical thickness is increased), all other param- 
X R 

eters remaining fixed, the value of q decreases and becomes asymptotic to zero. 

This trend is the same trend exhibited in the problem of radiant-energy transfer between 
two plane parallel walls separated by an absorbing and emitting, but motionless and non- 
heat -conducting gas. (See, for example, ref. 11.) Increasing y while the other param- 

R 

eters are held fixed results in an increase in q . This trend is the reverse of that for 
a transparent shock layer. (See fig. 21.) The results of the small-perturbation analysis 
(see fig. 9) show that this reversal occurs at intermediate values of the Bouguer number 
kp. Finally, it is apparent from figure 26 that an increase in surface reflectivity r w 
for fixed values of the other parameters results in a decrease in q w . The change in q 
with r w satisfies the inequality 



which agrees with the physical argument presented in the discussion of the small- 
perturbation solution. 


RADIATION -DEPLETED SHOCK LAYER 


The Strong Radiation Approximation 

When the radiation cooling parameter e is very much greater than both one and 

2 

kp, the Bouguer number squared, a particle leaving the shock with an initial enthalpy of 
1 2 

— W M will emit radiation at such a rapid rate that it will lose a significant amount of its 
energy before traveling the distance of a photon mean free path. Because this energy is 
emitted in a region of small optical thickness adjacent to the transparent shock, most of 
it will escape from the shock layer, and the enthalpy level within the shock layer will be 
very small in comparison with the value at the shock. In fact, as will be shown subse- 
quently, the zero-order solution for the enthalpy in the interior of the shock layer is 
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identically zero. It is for this reason that the term "radiation-depleted shock layer" 
has been coined. Of course, the narrow region adjacent to the shock in which the large 
change in enthalpy occurs can be described as a boundary layer and boundary-layer tech- 
niques can be applied to obtain solutions in it. 

The conditions which must hold in order that there be a radiation-depleted shock 
layer (that is, e very large and kp not too large^ occur at high altitudes for rather 
large objects (shock radius greater than 1 meter) entering at extremely high speeds 
(entry speeds close to 70 km/sec). It is not at all clear that the requirement for chemi- 
cal equilibrium can be satisfied under these conditions, particularly in view of the exis- 
tence of a shock boundary layer in which large changes occur over a short distance, and 
hence, a short time. Despite this objection, the solutions for the radiation-depleted shock 
layer represent an interesting limiting case which should lead to an increased under- 
standing of the radiating shock layer and provide a firm base for extension into areas of 
more practical concern. 


Analysis 


Once again, as was the case for the optically thin and optically thick shock layers, 
analysis can be facilitated through the use of the substitute kernel approximation. In 
this case, the energy equation, written in terms of the optical path length, is 

[f(r) h'(r)] M - | eB”(r) - | k|f(r) h’(r) = 0 (206) 


Here f is the nondimensional stream function, h the nondimensional enthalpy, B the 
nondimensional blackbody emissive power, e the radiation cooling parameter, kp the 
Bouguer number, and r the optical path length. It should be remembered that use of 
this equation restricts the analysis to gray gases only. The boundary conditions on the 
enthalpy are as before 


h(r A ) = 1 


(207) 


and the integral condition 


f(r) h’(r) - | e 


B(t) 
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B(t) e 


|k p |t-i 


dt 



P r w 


e 





(208) 


, \ 3 9 k P T 

The particular form of the substitute kernel employed here, that is, E2(kpTj ~ - e * 
was chosen for simplicity. Somewhat greater accuracy might be achieved by letting the 


62 


constants depend on the optical depth k p r A . However, it was not felt that this procedure 
would lead to a better understanding of the radiation-depleted shock layer. 

The momentum equation, in terms of the Dorodnitsyn coordinate, and the boundary 
conditions on the nondimensional stream function are 

2i{rj) f"(? 7 ) - [f’(7?)] 2 + a 2 h = 0 

f(0) = 0 

f W = 1 

When the radiation cooling parameter e is very much greater than one and very 
2 

much greater than k p , the energy equation (eq, (206)) admits the asymptotic solution 

B(t) = C 1 + C 2 t (213) 

Substitution into the asymptotic form of the integral condition (208) reveals that each 
of the two constants must be identically zero. Thus, in the interior of the shock layer 
B(t) and h(r) are zero. In this case, the density is infinite and the momentum equa- 
tion can be readily solved for f ( 77 ) with the result that 



and, of course, the shock standoff distance tends to zero. 

In order to investigate the shock boundary layer, it is convenient to introduce the 
"stretched" coordinates 

s = ( t a - ? y n (215) 

and 

? = (V A - V)e n (216) 

into the energy and momentum equations, respectively. Performing the usual manipula- 
tions (details are presented in appendix C) shows that the boundary layer is characterized 
by the parameter e * and it would seem proper to expand both the boundary-layer and 
asymptotic solutions in power series of this small parameter. However, f a (7?) (where 
the subscript a indicates the asymptotic solution valid far from the shock) is not analytic 


(209) 

( 210 ) 
( 211 ) 

( 212 ) 
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in e"* near e - ^ = 0, but is analytic in e~^^. Consequently, the solutions must be 
expanded as power series in e" 

The lowest order form of the energy equation in the boundary layer is 


^ B b,0( 9+ f B b,0 (a = C l +C 2« 


(217) 


The subscript b,0 has been used to denote the zero-order boundary-layer solutions. 
This equation must satisfy the boundary conditions that both B b g(|) and q(|) 
vanish as | increases without limit. Thus, the constants and C 2 are’ both 

identically zero. The third condition to be satisfied is 

B b0 (0) = 1 (218) 

The solution to equation (217) subject to the boundary conditions is 



(219) 


Solution of the momentum equation in the boundary layer gives the zero-order form 
of the nondimensional stream function 


f bj0 <£) = 1 (220) 

These zero-order solutions can be used to generate solutions of higher order. 
Mathematical details are presented in appendix C. In general, the analysis follows the 
procedure outlined by Van Dyke. (See ref. 38.) 

A complete listing of these solutions up to second order in the parameter e“ 
is presented. 

Zero-order solutions.- The zero-order solutions are: 


B a,0« = 0 



_ 2 C / dh\ dB 

3 Jr \dB/ B 
w b,0 


( 221 ) 

( 222 ) 


(223) 


f b,0^> = 1 


(224) 
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Vo = — 2X a ~ ^ = 1 + V 2 x(l - x) 


T Xo 

A ’° /Cp^O) 

First-order solutions .- The first-order solutions are: 

B a,l< T > “ 0 


B a,2 (0 > 

B b>1 ( £) = 0 
’ 0 


'A,0 


1 - 


Vo, 
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v> 

(226) 

(227) 

(228) 

(229) 

(230) 

(231) 
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f a,2 (??) = 
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KpW 

Second-order solutions .- The second-order solutions are: 
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A = 3 l ^j kpTa .° 


(235) 
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Radiant heat flux and standoff distance.- The total radiant heat flux to the wall q: 


R 

w 


and the ratio of the shock standoff distance to the shock standoff distance for radiation- 
less flow A are given by the following expressions: 

~ R 

(238) 
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Results and Discussion of Solution for Radiation-Depleted Shock Layers 

In the analysis of the preceding section, it was convenient to use the black-body 
emissive power B rather than the enthalpy h as the dependent variable. This choice 
necessitated the assumption that the thermodynamic and optical properties (in particular 
h and k^J be analytic functions of B in the interval (0,1). Unfortunately, this condi- 
tion does not hold for the correlations of the section "Stagnation Model for a Radiating 
Shock Layer" when written in terms of B in the limit as B approaches zero. 
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This difficulty can be circumvented through the use of analytic approximations to the cor- 
relating functions. For example, the enthalpy might be approximated by the function 


h=(B + B *) 1/ ' 2 (240) 

where B* is a very small positive number. Use of formula (240) in the expres- 
sion (223) results in the following solution for B^ q(|): 


where 


B b,0^ ~ 


-3\/b* V 
1 + Ce 3VB g j 

,1 - Ce- 3ji V 



c = \l 1 + B* - \/B* 
\l 1 + B* + \Jb* 


(241) 


(242) 


For large values of that is, 

k »|i/b* 

the value of B^ q(|) is directly proportional to B*. It is clear then that B* should 
be chosen sufficiently small to insure that B^ q(£) is nearly independent of B* for 

values of B^ q as small as e~^B a 2 ( 0 ). 

Because of the unlikelihood of establishing local thermodynamic and chemical equi- 
librium in a physical shock layer under those conditions for which this model analysis is 
supposed to apply, it would be somewhat superfluous to present the results of detailed 
calculations for the enthalpy distribution and shock standoff distance. It is sufficient to 
say that the shock layer is characterized by an enthalpy boundary layer immediately 
behind the shock the thickness of which (in terms of the Dorodnitsyn variable rj) is char- 
acterized by the inverse of the radiation cooling parameter e - ^. It should also be 
pointed out that the shock boundary layer is always very much less than a photon mean 
free path and hence is always optically thin. The enthalpy level in the interior of the 
shock layer is of order of magnitude kp^/e . The ratio of the shock standoff distance to 
the shock standoff distance for radiationless flow is of order of magnitude e - *. 

Curves representing the magnitude of the radiant heat flux which is absorbed by the 
wall q B are presented in figure 29. In the optically thin limit (kp-r^ q « l), the radiant 

heat flux approaches the "available energy limit" of (l - r w )y2. As the optical thickness 
of the wall increases and absorption becomes more important, less of the energy emitted 
from the gas in the shock boundary layer in the direction of the wall is able to penetrate 
the shock layer and reach the wall before being absorbed. Part of what is absorbed is 
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then reradiated in the forward direction and escapes from the shock layer through the 
transparent shock. Finally, as kpT^ q tends to infinity, none of the energy emitted in 
the shock boundary layer reaches the wall and the radiant flux incident on the wall 
vanishes. 


RADIATING SHOCK LAYERS 


Discussion of the Approximate Solutions 


Four different approximate stagnation-point solutions for an inviscid, radiating 
shock layer have been obtained. Each one represents a limiting case for some combina- 
tion of the radiation cooling parameter e and the Bouguer number kp. The regions of 
validity of the approximate solutions are depicted in figure 30. The boundaries as drawn 
pertain only to a gray gas with constant absorption coefficient. It would be necessary to 
redraw the boundaries for each nongray gas and for every change in the enthalpy depen- 
dence of the absorption coefficient. As was pointed out in the section "Optically Thin 
Shock Layers," the criterion for validity of the optically thin solution is that the gas layer 
be optically thin in all wavelength intervals in which a significant amount of energy is 
transported by radiation. For a gray gas, this criterion means kpT^ « 1. Thus, the 
boundary is not specified completely by kp but varies with e as well ^since kpT^, 
the shock-layer optical thickness depends on e as well as kp when the absorption 
coefficient is a function of the enthalpy). When applied to a 'nongray gas, the criterion for 
validity of the optically thin solution is always more restrictive than the condition that the 
Planck mean optical depth be small. Thus, the boundaries for all nongray gases are dis- 
placed to the left of the boundary for the "Planck-equivalent" gray gas. A "Planck- 
equivalent" gray gas is one in which the wavelength-independent absorption coefficient is 
equal to the Planck mean absorption coefficient in the nongray gas. 


The location of the boundary for the small-perturbation approximation depends on 
the radiation cooling parameter e and the enthalpy variation of the absorption coefficient. 
The value of e for which the solution will yield results of a given accuracy is reduced 
with an increase in y (where y is the exponent in the correlation formula Kp = h^), 
because of the reduced accuracy of the truncated expansion for Kp. Since the small- 
perturbation solution was shown to be correct to second order throughout most of the 
domain of the problem, the condition for establishing the boundary is « 1. (This 
condition holds when y = 0. When y = 4, the proper condition becomes lOe « l.j The 
location of the boundary does not depend on the wavelength dependence of the absorption 
coefficient. 


The most restrictive condition limiting the validity of the optically thick analysis 
for moderate values of the radiation cooling parameter e is the thickness of the enthalpy 
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layer adjacent to the shock. This thickness is characterized by the inverse of the Bouguer 
number l/kp. 

condition (e/kf 

boundary. This latter condition insures that the enthalpy boundary layer adjacent to the 
wall is thin compared with the shock standoff distance. The optically thick analysis is 
restricted to the case of a gray gas but could be extended rather simply to the case of 
an absorption coefficient with a step-function dependence on wavelength for which the 
step heights are either a(h) /Cp(h) or zero. There is no restriction to the number or 
width of the steps. The only changes that would appear in the formulas would be the 
substitution of a(h) /cp(h) for Kp(fc) and B(h)/o;(h) for B(h). The boundary to the 
region of validity of the optically thick shock-layer analysis would-be displaced to the 
left for this particular class of nongray gases. 

The region of validity of the radiation-depleted shock-layer analysis is restricted 

2 

by the conditions e » 1 and e » kp. The first condition insures that the thickness of 
the enthalpy boundary layer adjacent to the shock is small compared with the shock 
standoff distance, whereas the second condition insures that radiation is the preponderant 
mode of energy transport within the shock layer. Like the optically thick analysis, the 
radiation-depleted shock-layer analysis is restricted to gray gases but can be extended 
to the nongray model absorption coefficient with multiple steps of uniform height. Use 
of such a nongray model would cause a leftward shift in the boundaries to the region of 
validity. 

In order to relate the radiation shock-layer regimes to the problem of entry into 
the atmosphere of the earth, several trajectories are indicated on the e - kp map pre- 
sented in figure 31. The arrows indicate the direction of increasing time. Trajector- 
ies 1 and 2 represent iron spheres of radius 1 meter and 1 centimeter, respectively, 
entering vertically with an initial velocity of 70 km/sec: (No account has been taken of 
mass loss of these spheres due to ablation.) Trajectories for all other objects of the 
same size and lesser or equal initial velocities must lie below them in the e - kp space. 
The third and fourth trajectories correspond to the entry of round-trip Martian probes of 
different nose radius which would encounter some of the more severe heating conditions 
of the currently envisioned class of manmade objects. It is apparent that the small- 
perturbation approximation has considerable utility for the analysis of radiation effects 
on the entry of manmade objects. (See ref. 47.) It also appears that the optically thin 
shock-layer analysis might enjoy wide applicability. However, in the more realistic case 
of a nongray gas the boundary would be shifted to the left one or two orders of magnitude 
in the Bouguer number kp, and would thereby considerably reduce the practical useful- 
ness of the optically thin approximation. The optically thick and radiation-depleted shock- 
layer analysis would seem to be nearly devoid of direct practical usefulness, both because 
of the inaccessibility of the proper magnitudes of the parameters e and kp to objects 


Thus, the criterion for validity is kp » 1. For larger values of e, the 
\l/2 

« 1 becomes more restrictive and must be used to establish the 


69 


of interest and because of the restriction of these analyses to the gray case (and the 
simple nongray model absorption coefficient with multiple steps of uniform height) . 


A Model Earth-Entry Environment 


The four approximate solutions can be used to compute the radiant heat transfer to 
the stagnation point over a wide range of the radiation cooling parameter e and the 
Bouguer number kp. The results depend on the particular gas, the surface reflectivity, 
and the size of the object and must be recomputed for every change in these variables. 
Actually, the size of the object is important only if the exponent y (which appears in the 
correlation formula K-p = h^ varies throughout the e - kp space. In this event, the 
value of h at which a change in y occurs depends on the parameter kp = P g Kp g A^ 
which is influenced by the body size through the radiationless shock standoff distance A^. 

Contours of constant values of T s (the temperature immediately behind the shock) , 
P s k p (the Planck mean volume absorption coefficient immediately behind the shock), 

e /kp, and x (the ratio of free-stream density to the density immediately behind the 
shock) on plots of ambient density ratio P^/Pq as a function of free-stream velocity 
up to 70 km/sec are presented in figure 32 for a model earth-entry environment. 

This entry environment was obtained by combining the thermodynamic and optical prop- 
erty correlations presented earlier with the strong normal-shock relations. The resulting 
formulas are: 



p s*P,s = ( 7 * 94 X i 0 " 26 )^- 1 p^J T 
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for the lower temperatures (less than about 20 000° K) and 
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for the higher temperatures. 
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The values of x and y (the exponent in the correlation formula 
vary greatly over a rather extensive range of ambient densities and velot. 
quently, it was decided to fix these quantities at the constant values, y = t 
y = 4.0, for the discussions which follow. 



Radiant Heat Transfer 

The rate of radiant heat transfer to the stagnation point of a blunt object was 

calculated by the four approximate methods for a wide range of the radiation cooling 
parameter e and the Bouguer number kp. The results are presented in figure 33 as a 
plot of q w against kp for various values of the ratio e/kp. This ratio, sometimes 
known as the inverse of the Boltzmann number, was used because it is what might be 
termed an "environmental parameter," that is, a parameter dependent only on free-stream 
conditions (ambient density and velocity) and not on body geometry. The Bouguer number 
kp, on the other hand, is directly proportional to the body nose radius for a given set of 
free-stream conditions. Thus, each curve in figure 33 can be thought of as representing 
the effect of body nose radius on radiant heat transfer at a given trajectory point. 

For the purpose of calculating the results presented in figure 33, the shock-layer 
gas was assumed to have a gray mass absorption coefficient which varies as the fourth 
power of the enthalpy. The surface of the object was considered to be nonreflecting. 

The dashed portions of the curves do not represent computed data, but rather represent 
arbitrary connections across regions in which the various approximate solutions are 
invalid. 

The radiation cooling parameter e is equal to the radiant flux which leaves each 
side of a transparent, isenthalpic gas slab in which the nondimensional enthalpy takes the 
value one. Hence, this product represents an upper bound to the rate of radiant heat 

r> 

transfer to the stagnation point (or wall) q When this product is small, the rate of 
energy loss through radiation is small, and the average intensity is only slightly per- 
turbed from the isenthalpic value. However, as e increases (the Bouguer number kp 
remaining very much less than one), the increased energy lost by radiation results in 

O 

decreased levels of enthalpy and average intensity. Hence, q w becomes a decreasing 
fraction of e. Finally, as e becomes very large (kp still small), nearly all the energy 
is removed from the shock layer by radiation and q^ which represents the rate at which 
radiant energy leaves one side of the transparent layer approaches the physical maximum 
of 1/2. 

As kp increases toward and beyond unity, absorption becomes important and this 
mechanism, which tends to inhibit radiant energy transfer, halts the increasing trend of 
q w with kp. As kp continues to increase, the trend is reversed and q^; decreases 
and becomes asymptotic to zero. Consequently, the curves of rate of radiant heat transfer 
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to the stagnation point q“ against Bouguer number kp for constant values of the 
ratio e/kp have maximums the locations and heights of which depend on e/kp. It can 
be inferred from this statement that for every altitude and velocity in this simple model 
atmosphere, there is a finite value of nose radius for which the rate of radiant heat trans- 
fer to the stagnation point will be a maximum. It should be kept in mind that this analysis 
assumes a gray gas and that absorption in the free stream ahead of the shock has been 
neglected. 

In order to obtain some understanding of the effects of radiation cooling, gray 

absorption, and spectral absorption on the rate of radiant heat transfer to the stagnation 

point, a series of calculations utilizing various approximations were performed. The 

results of these calculations corresponding to a free-stream velocity of 14.2 km/sec and 

an altitude of 32.4 km are plotted against body nose radius R;n in meters in figure 34. 

The curve labeled 1 was computed by assuming that the shock layer was both isenthalpic 

and nonabsorbing. In this case the rate of radiant heat transfer to the stagnation point 

R 

is given by the simple expression q = e . This approximation was used in the early 
estimates of radiant heating. (See refs. 1 and 2.) Curve 2 was computed by assuming 
that the shock layer was isenthalpic and contained a gray, absorbing gas. The effect of 
gray absorption is seen to be small (under the conditions of this example) for a nose 
radius as large as 0.1 m. The third curve was obtained using the transparent approxi- 
mation discussed in the section "Optically Thin Shock Layers." This assumption of a 
nonabsorbing but radiation- cooled shock layer is frequently employed in the literature. 
(See, for example, refs. 3 to 7.) For this example, at least, the effect of radiation cooling 
is more important than the effect of gray absorption for nose radii of 0.1 m or less. 

Curve 4 contains the effects of both radiation cooling and gray absorption. These com- 
bined effects are included in the numerical solutions of Howe and Veigas (ref. 9). It can 
be seen that for small nose radii (less than about 0.1 m), gray absorption has little effect. 
However, gray absorption plays an increasingly important role as the radius increases. 

The final curve, curve 5, includes the combined effects of radiation cooling and nongray 
absorption. The absorption coefficient used in these calculations was the step function 
model introduced in the section on the small-perturbation solution. (See fig. 13.) The 
curve is limited to small values of nose radius because of the restricted region of valid- 
ity of the small-perturbation method with which this curve was computed. It is very 
apparent from these results that nongray effects cannot be ignored if one wishes to obtain 
a realistic evaluation of the radiant heating of objects during entry at hyperbolic velocities. 

The analysis of this paper has been restricted to a shock layer with plane-parallel 

geometry. The largest effect of this assumption is felt in the calculation of the rate of 

radiant heat transfer. Koh (ref. 21) has shown that the plane -parallel geometry assump- 

R 

tion can lead to an overestimation of q w by no more than 15 percent when the gas is 
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transparent to its own radiation and when the shock standoff distance to shock radius 
ratio is no greater than 0.05. As the Bouguer number kp increases, the size of the 
error decreases and vanishes when the shock layer becomes optically thick. Because 
the effective optical thickness of a nongray shock layer is greater than that for a Planck 
equivalent gray gas, the error due to geometry will be smaller for a given Bouguer num- 
ber in the more realistic nongray case. 

Convective Heat Transfer 

Even though the analysis of this investigation is based on the assumption that the 
gas in the shock layer is inviscid and nonheat conducting, it is possible to draw some con- 
clusions regarding the coupling between radiant heat transfer and convection heating. The 
convective heating rate (sometimes referred to as the aerodynamic heating rate) is, in 
the case of a laminar boundary layer, the rate at which heat energy is transferred to the 
body surface by means of conduction. 

To first order in the boundary-layer parameter (Npe) - *^, the convective heating 
rate is proportional to the enthalpy difference across the conduction boundary layer. If 
the wall is cold (as has been assumed throughout this investigation), the enthalpy of the 
wall can be neglected and the convective heating rate becomes proportional to the enthalpy 
at the outer edge of the boundary layer. The location of the outer edge depends upon the 
Peclet number. Since it has been assumed throughout this investigation that the viscous 
boundary layer is thin (in terms of both the Dorodnitsyn coordinate and the optical path 
length), the location of the edge of the viscous boundary layer is arbitrarily specified as 
V/V^ = 0-05 for both the small-perturbation and the transparent solutions. The rapid 
change in enthalpy near the wall, particularly for the transparent approximation which 
gives a value of zero for the enthalpy of the gas adjacent to the wall, necessitates choosing 
an edge location other than zero. For the optically thick and radiation-depleted shock 
layers, it is more convenient to specify the edge of the viscous boundary layer in terms 
of the normalized optical path length r. The variation of enthalpy near the wall is very 
small in the case of the radiation-depleted shock layer. Consequently, the edge of the 
viscous boundary layer can be considered to be located at t = 0 for this case. A wall 
boundary layer due to radiation has been shown to exist in the optically thick shock layer. 
This wall boundary layer is always thicker than a photon mean free path and, of course, is 
very much thicker than the optically thin viscous boundary layer. Therefore, r = 0 can 
be considered as the edge of the viscous boundary layer for this case also. Values of the 
enthalpy h e at the outer edge of the viscous boundary layer have been determined from 
the four approximate solutions for a wide range of the ratio of the radiation cooling param- 
eter to the Bouguer number e/kp and the Bouguer number kp. The results are shown 
in figure 35. The dashed portions of the curves represent arbitrary connections across 
regions of nonvalidity. 
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The quantity h e is a rough approximation to the ratio of the convective heating 
rate for a radiating shock layer to that for a nonradiating shock layer. When radiant 
energy transport is important, the convective heating is reduced from the radiationless 
value (h e = l). The effect becomes larger as both e/kp and kp increase. It is 
interesting to note that the convective heating continues to decrease for increasing kp 
even when the shock layer is optically thick and the rate of radiant heat transfer is 
decreasing as a result of absorption. Even though the total heating rate (radiant plus 
convective) cannot be deduced from an inviscid analysis, it is apparent that the total 
heating rate decreases with increasing shock-layer optical thickness for all values of 
kp at least as large as the value for maximum rate of radiant heat transfer to the stag- 

D 

nation point q w . 

Of course, the results of figure 35 only give an order-of-magnitude estimate of the 
radiation- convection heating coupling. Not included are the effects of variable transport 
properties, enthalpy gradient at the edge of the boundary layer, and differences in the 
characteristic Reynolds and Prandtl numbers between the radiating and nonradiating 
cases. Also no account has been taken of the effect of radiation in the boundary layer. 

In the cooled region of the boundary layer adjacent to the wall, the gas will absorb more 
radiant energy than it will emit. This condition will tend to increase the slope of the 
enthalpy distribution adjacent to the wall and thereby increase somewhat the convective 
heat transfer. 

The effects of radiation cooling, gray absorption, and spectral absorption on the 
ratio of convective heating rate for a radiating shock layer to that for a nonradiating shock 
layer h e is shown in figure 36. It is apparent that radiation cooling plays the major 
role whereas absorption (both gray and nongray) tends to reduce the effectiveness of 
radiation cooling. The calculations for curves 1 and 2 ignored cooling. Consequently, 
no reduction in the calculated value of the convective heating rate was obtained. Curve 3 
includes radiation cooling and ignores absorption. Thus the reduction in the calculated 
value of the convective heating rate is maximized in this approximation. Finally, curves 4 
and 5 indicate that absorption inhibits the effectiveness of radiation cooling, and since 
absorption is more important in a nongray gas than it is in a Planck equivalent gray gas, 
the rate of convective heating will be greater in the nongray case. 

The Role of the Radiation Cooling Parameter and the Bouguer Number 

The radiation cooling parameter e admits of several physical interpretations 
which are useful in the understanding of the radiating shock layer. Of these interpreta- 
tions, one of the most useful is the following: 
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e = 


( Time required by element of volume\ 
to traverse distance A A at rate of I 

emergence from shock / 

2(Energy of element of volume upon emergence from shock) 


Here A^ is the shock standoff distance in a nonradiating (or adiabatic) flow. 

It can be seen from this interpretation that the radiation cooling parameter is indic- 
ative of the slope of the enthalpy distribution immediately behind the shock. In fact, in 
the transparent limit there is a direct relation between e and the initial slope; that is, 


(See, for example, fig. 17.) 



= 2e 


In the case of an optically thick shock layer, the initial enthalpy gradient is reduced 
by absorption. However, a lower bound to the gradient is the value e (one-half the 
transparent value) because the emergent elementary volume will emit at least twice as 
much energy as it absorbs; it emits energy at a rate proportional to the Planck function 
at the equilibrium shock temperature T s in both the upstream and downstream direc- 
tions while it absorbs energy at a rate at most proportional (by the same factor; the 
monochromatic volume absorption coefficient) to the Planck function at temperature T s 
from only the downstream side. 


A physical interpretation of the Bouguer number is: 


k p = 


Radiationless shock standoff distance, A^ 

Planck average photon mean free path in gas emerging from shock 


Only when conditions do not vary greatly across the shock layer will the Bouguer number 
be indicative of the Planck mean optical thickness and only when the gas is nearly gray 
will the Planck mean optical thickness be indicative of the various important monochro- 
matic optical thicknesses. Consequently, critical values of the Bouguer number are sub- 
ject to a number of influences; among them, the enthalpy and spectral variation of the 
absorption coefficient and the value of the radiation cooling parameter. For example, the 
value of kp for which absorption first becomes important is about 0.1, when the radia- 
tion cooling parameter is very much less than one and the absorption coefficient is inde- 
pendent of wavelength. When the absorption coefficient varies spectrally as shown in 
figure 13, and when e is very small, absorption begins to become important for Bouguer 
numbers as small as 0.001. With e about 10 for a gray gas, absorption is important 
for values of the Bouguer number greater than about three. Despite these drawbacks, the 
Bouguer number as defined in this investigation is about the best "a priori" indicator of 
the importance of absorption that can be obtained. 
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When the radiation cooling parameter e is very much less than one, an elementary 
volume of gas will lose very little of its energy by radiant emission in the time required 
to traverse most of the shock layer. (Of course, it takes an elementary volume of gas 
traveling along the stagnation streamline an infinite time to reach the wall.) Hence, 
radiation cooling of the shock layer will be slight. When the radiation cooling parameter 
is very much greater than one, an elementary volume of gas will emit energy at such a 
rapid rate that the energy of the volume will be reduced a significant amount before it 
leaves the vicinity of the shock. This is true whether the shock layer is optically thick 
or optically thin (that is, regardless of the size of the Bouguer number). This physical 
argument is used to establish the existence of the thermal boundary layer behind the 
shock in the radiation-depleted shock layer. If the shock layer is optically thick, the 
reduction in enthalpy will continue only so long as the elementary volume is within about 
a photon mean free path of the shock. Beyond this point, the elementary gas volume 
receives radiation from all sides and begins to establish a condition of radiative equi- 
librium with its surroundings. The energy lost during the time required by the ele- 
mentary volume to travel a single photon mean free path is characterized by the ratio of 
the radiation cooling parameter to the Bouguer number e/kp. (it was shown earlier that 
the enthalpy level in the interior of an optically thick shock layer was characterized 
solely by the parameter e/kp.) 

Within the interior of an optically thick shock layer, radiation heat transfer can be 

treated in a manner analogous to conductive heat transfer. Thus, one would expect that 

a parameter analogous to the Peclet number could be constructed which would suggest the 

nature of the enthalpy boundary layer adjacent to the wall. Such a parameter, which is a 

ratio representing the relative importance of convective to radiative heat transfer is 

2/9 

given by the grouping kp/e. Since the thickness of the enthalpy boundary layer is char- 
acterized by (Np e ) - ''' / ^ in the conduction problem, one expects, by analogy, the thickness 
of the enthalpy boundary layer adjacent to the wall in an optically thick radiating shock 
layer to be characterized by £ V2/ kp . 

The importance of this parameter (in a somewhat 
different form) and its analogy with the Peclet number was pointed out previously by 
Goulard. (See ref. 23.) 

The importance of the surface reflectivity r w depends on the importance of 
absorption in the shock layer. When absorption is negligible, the effects of surface reflec- 
tivity are negligible because all photons originating within the shock layer will escape the 
layer and it matters not whether some of these photons are absorbed by the cold wall or 


^In the optically thick shock-layer analysis, the gas was assumed to be gray. Thus, 
the fact that the Bouguer number was based on a Planck mean absorption coefficient was 
of no consequence. However, in the case of a nongray gas, it would probably be more cor- 
rect to use a Bouguer number based on a Rosseland mean absorption coefficient. 
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reflected by the wall into the free stream. However, when absorption is important, the 
reflected photons have a large probability of being recaptured in the shock layer. Thus, 
an increase in surface reflectivity tends to raise the enthalpy level of an absorbing gas 
in the vicinity of the wall. 

In this section it was shown that both the radiation cooling parameter e and the 
Bouguer number kp play prominent and interrelated roles in determining the character 
of the radiating shock layer. Further, it was shown that the spectral variation of the 
absorption coefficient greatly influences the role of the Bouguer number. In general, 
then, one cannot ignore either of the processes of radiation cooling and nongray absorption. 

CONCLUDING REMARKS 

A mathematical model for the stagnation region of a radiating shock layer was 
derived in this investigation subject to the following conditions: (1) the gas in the shock 
layer is in local thermodynamic and chemical equilibrium, (2) the body geometry is axi- 
symmetric, (3) there is no mass addition to the flow from the body surface, (4) the thick- 
nesses of the shock and the viscous boundary layer are small in comparison with the shock 
standoff distance, and (5) absorption in the free stream ahead of the body is negligible. 

The divergence of the radiant flux vector, which appears in the energy equation, was for- 
mulated to include a wavelength-varying absorption coefficient. The body surface was 
considered to be cold and to reflect diffusely and independently of wavelength a fraction 
equal to the surface reflectivity r w of the incident radiation. The results of a boundary- 
layer analysis indicate that the equations for the flow in the inviscid region are indepen- 
dent of the boundary-layer equations only when the boundary layer is optically thin or 
optically thick. It has been assumed throughout this study that the boundary layer is 
optically thin. Simple correlation formulas for the thermodynamic and optical proper- 
ties of high -temperature equilibrium air were developed and used herein. 

The general form of the governing system of equations was found to be integro- 
differential in character. The solution of this system is extremely difficult to find even 
with numerical techniques and high-speed electronic computing machines. The approach 
of this investigation was to take advantage of the simplified form to which the governing 
equations were reduced when the radiation cooling parameter e and the Bouguer num- 
ber kp took on limiting values and obtain approximate analytic solutions if available. 

It was found that the general problem reduced to a singular perturbation problem in each 
of the four cases studied. A small-perturbation solution valid when the energy lost to the 
shock layer by radiation is small (that is, when the radiation cooling parameter is small) 
was derived. The Poincare-Lighthill-Kuo (P-L-K) perturbation of coordinate method 
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was used to obtain a uniformly valid solution. This solution was used to study radiation 
cooling, absorption, effects of surface reflectivity, and effects of nongray optical 
properties. 

An optically thin shock-layer method of solution utilizes an expansion in terms of 
the Bouguer number kp to reduce the governing system to purely differential form. 
Again it was necessary to resort to the P-L-K method to obtain a uniformly valid solu- 
tion. This solution was used to study radiation cooling, absorption, and the effects of 
surface reflectivity. 

The optically thick approximation, valid when the optical thickness of the shock 
layer is very large (that is, the Bouguer number very much greater than 1) was used. 

The governing equations were reduced to differential form through the use of a substitute 
kernel approximation. Two thermal boundary layers were seen to exist; one adjacent to 
the shock and the other adjacent to the wall. It was noted that the Rosseland approxima- 
tion together with a properly specified temperature-jump or slip condition at the wall 
reduces the governing equations to the same form as the substitute kernel approximation 
in the interior or isenthalpic portion of the shock layer and in the wall boundary layer. 
However, the Rosseland approximation with slip conditions was found to be inadequate for 
analyzing the shock boundary layer. The optically thick solutions were restricted to gray 
gases but were used to study radiation cooling, absorption, and the effects of surface 
reflectivity. 

A radiation-depleted shock layer was analyzed. This approximation is valid when 
the rate at which energy is radiated away from the shock layer is nearly equal to the rate 
at which energy enters the shock layer (that is, the radiation cooling parameter is very 
large) so that the enthalpy level is very much less than the radiationless value. The sub- 
stitute kernel approximation was used to reduce the governing system of equations to dif- 
ferential form. The method of matching of inner and outer expansions was used to obtain 
solutions valid in the thermal boundary layer adjacent to the shock and in the interior of 
the shock layer. These solutions were restricted to gray gases but were used to study 
radiation cooling, absorption, and the effects of surface reflectivity. 

It is apparent from the results that radiation cooling first becomes important when 
the rate of energy lost by radiation from the shock layer is only about 1 percent of the 
rate with which energy enters the shock layer. Absorption in a gray gas begins to 
become important for shock layer optical thicknesses greater than about one-tenth. An 
increase in the surface reflectivity r w from zero reduces the radiant heat transfer by 
a factor of roughly 1 - r w and increases the heat -transfer rate to the wall by conduc- 
tion because of an increase in enthalpy level near the wall. 

The results of some nongray calculations obtained with the small-perturbation solu- 
tion are presented. The Planck mean absorption coefficient can be used to compute the 
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enthalpy distribution and the radiation heat-transfer rate to the wall as long as the optical 
depth of the shock layer is very much less than one in all wavelength regions in which a 
significant amount of radiant energy is emitted. For larger optical thicknesses nongray 
effects are very important. 

The various approximate solutions were used to compute the rate of radiant heat 
transfer to the stagnation point of blunt objects traversing an optically gray model earth 
atmosphere. The results of this computation indicate that at every altitude and velocity, 
there is a finite value of body nose radius for which the rate of radiant heat transfer to 
the stagnation point is a maximum. (This result is contrary to the earlier results, based 
on the assumptions of an isenthalpic and transparent shock layer, which indicated that the 
heating rate was directly proportional to nose radius.) A significant reduction in the com- 
puted value of the radiant heating resulted upon taking the nongray character of air into 
account. This result served to emphasize that the nongray character of gases plays a 
very real and important part in problems of radiation gas dynamics. 

In general, the coupling between radiant and convective heat transfer is such that 
increases in the rate of radiant heat transfer result in decreases in the rate of convective 
heat transfer to the body surface. Of course, the amount by which the total heating rate 
is affected cannot be determined from this inviscid analysis. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., February 18, 1970. 
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APPENDIX A 


THE VISCOUS BOUNDARY LAYER 


In this appendix, a boundary-layer analysis will be performed on the integro- 
differential system (52) to (59) to determine the form of the equations in the inviscid 
region and the viscous boundary layer and to determine under what conditions such a 
separation can be achieved. For convenience, the system will be rewritten here 


where 


2 \* 


i(r)) h’( 7 ?) + A 2 
Np r ? 2 (h) f”0?) 


|*5j(h) h’ (77) + el [77] = 0 

(Al) 

+ 21(77) f"0?) - [|(??)) 2 + a 2 h ( r ?) = 0 

(A2) 

f(0) = 0 

(A3) 

f'(0) = 0 

(A4) 

f K) = 1 

(A5) 

\ " a 
V ||2 X (1 - X) 

(A 6) 

h(0) = h w 

(A7) 

h ( V ) = 1 

(A8) 


4ot; 


€ = PcKp 


s^P.s A iT 


(A9) 


and A, introduced here for convenience of notation, is the inverse square root of the 
Peclet number. 

When the parameter A is very much less than one, a perturbation type solution 
can be attempted. However, the energy and momentum equations each lose the most 
highly differentiated term as A vanishes. As a result, neither the zero-order (in the 
small parameter A) solution for f ( 77 ) nor that for h(? 7 ) can satisfy all the boundary 
conditions. In particular, the conditions f'(0) = 0 and h(0) = 0 must be relaxed, and 


80 



APPENDIX A - Continued 


the perturbation solution will not be valid as 77 approaches zero. Thus, this problem 
is a singular perturbation problem of the "boundary layer" type. (See refs. 38 and 54.) 

In order to obtain the boundary-layer form of the equations, the "stretched" coor- 
dinate £ = X _a 77 is introduced where a is an as yet undetermined constant. It is also 
convenient, to avoid confusion, to introduce the change in notation 

i(£) = h(rj) (A10) 

g’(l) = f'fa) (All) 

J[S] = l[v] (A12) 

Equation (All) is written in this particular form because it is f'(rj) and not f(r]) which 
fails to satisfy the boundary condition at 77 = 0. 

When the stretched coordinate £, and the definitions (A10) to (A12) are introduced 
into system of equations (Al) to (A8), the only choice for a which will retain the most 
highly differentiated terms without loss of the most significant terms in the "unstretched" 
problem is a = 1. Thus, A and not A^ is the significant small parameter and the 
stretched coordinate is 

I = * -1 77 (A13) 

Perturbation solutions are now sought in the forms: 

OO 

i(£,A) = £ X\(J) (A14) 

n=0 

00 

g(«,X) = £ A n g n (ij) (A 15) 

n=0 

in the boundary layer, and 

OO 

h(77,X) = £ X n h n (7?) (A 16) 

n=0 

OO 

f(77,X) = ^ X%(7,) (A 17) 

n=0 

in the inviscid region. 
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APPENDIX A - Continued 


It shall be assumed that all functions of h (and i) are analytic about the value hg 
(and ig) so that they may be expanded in Taylor series about h = hg and i = ig in the 
following manner: 


F(h) = F^hg + Xh]^ + X 2 h2 + ..!)= F^hg) + XF^hgjhi + X 2 F^hg)h2 + ^ F|hg) h 


Dns 

oo 

(A18) 

i[n,x] = J 

n=0 

OO 

(A19) 

j[ 5 ,x] - 2 
n=0 

OO 
r — \ 

] (A20) 

a = ^- na n 

n=0 

OO 

(A21) 

= I X %,n 
n=0 

(A22) 


is also assumed without, for the present, specifying details of the terms Inj^7j and 

J nW- 

Furthermore, to insure compatibility of the boundary layer and inviscid solutions, 
it is necessary that the inner boundary condition on the outer solution be written in the 
form 

f(6) = 0 (A23) 

where 6 (the displacement distance) is specified by the matching condition 

lim g(£) = X -1 f(?7) (A24) 

£-*.00 

The quantity 6 depends on X and must be written in expanded form 
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OO 

6 = ^ X n 6 n (A25) 

n=l 

The term 6g was chosen to be zero because 6 is of order X. 

The system which describes the solutions valid in the inviscid region can be 
obtained by substituting expansions (A16) to (A22), and (A25) into system equations (Al) 
to (A8). The result is an infinite power series in X the sum of which is zero for all 
values of X. The only such series is one for which the coefficient of each of the X n 
terms is identically zero. These coefficients yield a set of recursive integro-differential 
systems. The system of zero order is 


*ofa) h ofa) + = 0 

(A26) 

2f 0 fo) fo’(r?) - 

, -| 2 o 

f 0 fa) + a^ofa) = 0 

(A27) 

o 

ii 

3 

"o 

SH 

(A28) 

f o( 

>.o) = 1 

(A29) 

' ' ^X(I-X) 

(A30) 

h O^A,o) = 1 

(A31) 


The system which describes the solutions valid in the boundary layer can be 
obtained by substituting expansions (A14), (A15), and (A18) to (A22) into system (Al) to 
(A8). As for the inviscid case, this procedure results in a set of differential systems. 
The zero-order system is 


j^o) i o(^) 


K(‘o 



+ g 0 (S) ij(«) + eJ 0 [«}» 0 
|2 


2g 0 (£) g 0 (!) - 
g 0 (0) = 0 

gj>(0) = 0 


Bn® 


+ a oi 0 (S) _ 0 


(A32) 

(A3 3) 

(A3 4) 
(A3 5) 
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lim gUt) = f^O) 
1—00 U U 


(A3 6) 


i 0 (0) = h w 


(A3 7) 


lim io(^) = h 0 (°) (A38) 

4-oo 

The first-order term of the displacement distance is found from the matching 
condition 


lim g 0 (4) = ^ - S^f^O) (A39) 

It is apparent that the zero-order solution for the boundary-layer equations depends 
only on the inviscid enthalpy level in the vicinity of the wall and not the enthalpy gradient. 
The enthalpy gradient will, of course, have an effect on the first-order boundary-layer 
solution. Thus, if the enthalpy gradient is very large, as it can be for a radiating shock 
layer, the boundary-layer solutions must be carried out to first order in X. When 
absorption is neglected, the enthalpy gradient at the body surface can be infinite and the 
analysis presented herein will not be valid. Burggraf (ref. 55) has shown how to over- 
come this difficulty. 

The divergence of the radiation flux (see eq. (64)) includes integrals which extend 
over the whole domain of the problem. It is convenient to separate each of these inte- 
grals into two integrals as follows: 


C T X A / \ / 


\ r T A,A 

/ \ 

J 0 B x (t,)E l( k p 

t a - ^ 

) dtk = J T * Bk 

K^J 


E l( k P T x " t'll dt 


and 



E 1 kp t a 



ds A 


(A40) 



E 2^ k P t x) 


dt-x 


l x 




where 



(A41) 

(A42) 
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is the monochromatic optical thickness of the boundary layer and £* is the thickness of 
the boundary layer in terms of the stretched Dorodnitsyn variable ij. 

It is convenient to redefine the monochromatic optical path length as 


t a = k a&> d £ 

VII 

(A43) 

t a = dr l + t a 

( T * > ^ 

(A44) 


In order to expand equation (64) as a power series in A, it is necessary to expand 
the exponential integral functions and all functions of h (and i) as well. Expanding the 
optical thickness yields, for 


t A = Ao A 



d £ + ^ il(£) «x(iq 


d£ + 


(A45) 


and for 7\ > t* 



The exponential integral functions can be expanded in the Taylor series 


(A46) 


E n (x - y) = EJ x 0 - y 0 + M*i - y i) + 


s E n (x 0 - y q) - - yi)E n _i^x 0 - y 0 ) + . . 

(A47) 


If the argument is of order A, 


E n (Ax) = E n (0) - A lX E n _!(0) + . . . + (-l) n - 2 A n “ 2 E 2 (0) 


(n - 2): 


(-l) n_1 A n_1 log e A 


x 


n-1 


(n - 1)! 


(-l) n-1 A n-1 x 


n-1 


(n - 1)1 


y + log e x 


+ . . . (A48) 


where y is Euler's constant (y = 0.577216). Use of this expansion, while it avoids any 

dependence of the terms In[^] on introduces terms of order A log e A into the 
boundary-layer solutions. 
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Incorporating the various expansions into equation (64) and separating the result 
into powers of X and X log e X yields the zero-order expressions: 



-2/fp[ho(r?)j Bjh 0 (?])] + k P k x[M^)]|JJ A K x[ h ofa T )] B x[ h 0^')] E l^ k p| T xfa) 
- T x( ? 7 , )|j d77» + 2 r w E 2 (kpT x (r 7 )) K x [h 0 (r7’)] B x [h 0 (r?’)] ^(kp^’)) drX dX 


(A49) 


and 


JoW = -2/<p[ i o(l)] B[i 0 (|)] + k p « x [ i o^)]|jJ A K x[ h oM B x[h 0 ^’)] E^kpr^')) drj’ 
+ 2r wj' T?A K x[ h 0^')J BxfhQ^’)] E 2 (kpT X (v')j drX dX (A50) 


The second of these expressions, which contains only definite integrals, is valid 
only when expansion (A47) holds. But equation (A47) converges in the first few terms 
only if the argument Xx ^or in the terms of this problem Xkpoj is small compared 
with 1. Thus, expressions (A48) and (A49) can be used only when the boundary layer is 
optically thin; that is. 


Xkpo* = kpT* « 1.0 


(A51) 
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MATHEMATICAL DEVELOPMENT FOR OPTICALLY THIN SHOCK LAYERS 


With the substitute kernel approximation, the divergence of the radiant flux vector 
can be written 

r°° CrV a - 2 k T >|T 1 ( 7 j)-T > (^)[ 

= l[r?] = 2 k p (? 7 ) B(t?) - 2kpj o K xfa)< j fl A k x {£) B x (£) e 1 d| 

1 ♦ r w e - ^ B x ( 9 e- 2 k P^ ( « d { ) dx 


(Bl) 

The monochromatic optical path length kpT x (rj) is given by the expression 


kp T x (7?) = k p yj K X U) d£ 


(B2) 

where k p is the Bouguer number 



k p = P s k P,s A A 


<B3) 

The approximate governing system presented in the section 
Layers" is 

"Optically Thin Shock 

i(v) h’(?7) + el[f?] = 0 


(B4) 

2 f ( 77 ) f "( 77 ) - [f'( 7 )j + a^h = 0 


(B5) 

f( 0 ) = 0 


(B 6 ) 

r-l 

II 

^<P 


(B7) 

f’(r? A ) - 

V ' 1 + ^2 X (1 - X) 


(B 8 ) 

h M = 1 


(B9) 
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The quantity ii is defined by the expression: 

h = — T A h(«) (BIO) 

J 0 

The Conventional Perturbation Procedure 

If the functions h/? 7 ;kp) and f^;kpj are analytic in the vicinity of kp = 0 (it 
is assumed herein that such is the case), they may be written in the expanded form: 

OO 

h (n;k p ) = ^ kph n (?7) 
n=0 

CO 

f(??;kp) = ^ kph n (r7) 
n=0 

It is anticipated that the first few terms of these expansions will provide an accurate 
estimate of the solution to the system equations (B4) to (BIO) when kp is small com- 
pared with unity. 

In addition, all quantities which depend on the parameter kp either directly or 
indirectly must be expanded in terms of kp. For example, a function 2^h(r?)j becomes 

jt[h(? 7 )] = “3 [h 0 (h) + kphjO?) + . . .] =^jh 0 (?7)] + k p 3^[h 0 (?7)J hj(r?) + . . . (B13) 

jf[h(77)] = 9 ^( 77 ) + kp'3 £ 0 (?7) hj(?7) + . . . (B14) 

The quantities 2q(v) and 3 ^( 77 ) have been introduced to simplify the notation. The con- 
stant 77 ^ is given by the expansion: 

= ^A,0 + k P^A,l + • • • (B15) 

Substituting the expansions (Bll) to (B15) into system equations (B4) to (BIO) yields 

f o 07) hg (??) - 2eK p , 0 (77) B 0 (? 7 ) + k p jfgfa) h ^ ( 77 ) + fjfa) h' 0 (v) - 2e K p ^(v) ^q(v) 

+ Kp jO 0?) B 0 (v) hi(??) + 2e(l + r w )J K X,0 ^ B x,0^ dxj + . . . = 0 

° " ° J " (B16) 


(Bll) 

(B12) 
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where 


]2i 0 (v) io(v) - + a 2 hoJ + k p j^fgfa) i^q) - 2i' Q (q) f ^(77) 

+ 2ig(q) f 1 (?7) + a 2 hij + . . . = 0 


f 0 (0) + kpfj(O) + . . . = 0 


f o(%o) + k P f l(^A,o) + ^A,l f o(^A,o) 

*0(^,0) + k p[ f l(^A,o) + ^A,l f o(^A,o)J + * • ‘ ■ 

^(^,0) + k P h l(^A,o) + ^A,l h o(^A,o) 


+ • • • — 1 


1 + \| 2 x(l - x) 


+ • • • — 1 


hQ + kphj + 


• = ^ a ’ 0h ° <ftd{+kp (^< o a, ° hiK)di 

- ^ C A, ° i> 0 <{) »( * + • • • 

<0 0 A ’°j 


Since the small parameter is arbitrary, equation system (B16) to (B21) can be 
satisfied only if each coefficient of each term is identically zero. This condition leads 
to a recursive set of purely differential systems. The zero-order system is 


h 'o(v) ~ 2 eK Pj0 (?]) B q (?7) = 0 


2f o 0?) fo'O?) - fgO?) 2 + a 2 h 0 = 0 


f 0 (0) = 0 


f o(^A,o) = 1 

1 + J2 X (1 - X) 
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H)(^a,o) 


= 1 


e 0 =— r a '° h o<«> d « 

° ”4,0 J 0 


The solutions to this system are 


C 1 dh _ 2e7? A,0 (1 - a*)x 

Jhn *p00 " * 

f 0 (7j) = /l - a\^r7 Aj0 


+ ar 


i 0 K p (h) B(h) a * ' ° e x 

„ \2 


\>,oy 


P^A, 




f 0 (??) = (1 - a*)x 2 + a*x 


^A,0 


l+\/2x(l - X ) 

1 + ^2h Q x(l - x) 


: o = y o ‘»0« 


dx 


where the definitions 


= a \|%?A,0 


x = 


V 


VA,0 


have been added to simplify the notation. 
The first-order system is 


f Q (v) h’i(??) - 2e k P j0 (??) B 0 (t7) + /c Pj0 (r?) B 0 (rj)Jh 1 (tj) 

= -tjW tyn) - 2e(l + r w )y” K xfi {v)^ o A ’° K Aj0 (4) B Xj0 (9 d{ 


d\ 


(B28) 

(B29) 

(B30) 
(B3 la) 
(B3 lb) 
(B32) 

(B33) 

(B34) 

(B35) 

(B36) 
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toiv) f-iiw) - f()(y) t'lln) + fo'C 7 ?) f i^) + | a2 ^i = 0 ( B37 ) 


fi(0) = 0 


(B38) 

f l( ?7 A,o) = ^a/o(^A,o) = ~ 

2b 4,1 

(B39) 

1 + \[2 X (1 - X) 
2(1 - a*)77 Ajl 

tl^Afl) = ■^A.l^A.o) = " 

(B40) 

”1,0 

h l(^A,o) = “^A.lbQ^A,^ 

1 = - 2e "A,l 

(B41) 


The solutions to this system are 


hi(7?) = -2e/c p ^ 0 (x) B 0 (x)<^?7 A x 


^A,O h l 


* a*2 N 

1 - a* + a 


(1 -x) 


2e ( 1 + r w) K P,0 (x ) B oWvl >0 § 


% ' (1 - a*)x + a* a 

poo r i 

i J„ *x,o«)j„ %o«'> B », <,(«’> «' 




dX 


f 0 (?) «p o (?) B o(?) 


d| 


(B42) 


fj(x) = x(2 - a*x) 


?7 


A,0 


7 a,i 


a*h. 


4h r 


7 A,0 


h l 



(B43) 


(B44) 


(B45) 


The expression for hj(?7) (eq. (B42)) can be simplified somewhat by transforming 
the integral over £ into an integral over hg through the use of equation (B23). The 
result of this transformation is 
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’J A,0 h 11 


h x (x) = -2e/c Pj0 (x) B 0 (x) \ v? A;1 + 


0 /L 


a*2) (1 - x) i (1 - a*)x + a* 

1 - a + — ;r —I ±r 10g„ -5 

2 / (1 - a*)x + a* a* e x 


(l + r w ) Kp;0 (x) B q (x) K xfi (V) B X>0 U') d«' ^ - 


1 K*(h) dh 


dA 


[/c p (h) B(h)]‘ 


(B46) 


The Poincare-Lighthill-Kuo Solution 

It can be seen on careful inspection of equation (B46) that the first-order term 
hj(x) displays a singular behavior in the vicinity of the wall (x = 0). Consequently, the 
assumed expansion for h(x;kp) diverges as the origin is approached and the perturba- 
tion solution is not uniformly valid. However, if the coordinate x is perturbed, the 
solution can be made uniformly valid. Thus, according to the P-L-K method, 

x = y + kpXj(y) + . . . (B47) 

where y is the transformed variable. The enthalpy when expanded in terms of kp 
with coefficients as functions of y, not x, becomes 

h(x;kp) = hj(y) + kph^y) + . . . (B48) 

and the nondimensional stream function 

f(x;kp) = fj(y) + kpf J(y) + . . . (B49) 

According to Pritulo (ref. 49), the coefficients in the P-L-K expansions can be related to 
the coefficients in regular expansions in the following manner: 


hj(y) = h Q (y) 

(B50) 

hj(y) = h 1 (y) + xj(y) h Q (y) 

(B51) 

f$(y) = f 0 (y) 

(B52) 

fi(y) = fj(y) + xj(y) fg(y) 

(B53) 
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The arbitrary quantity Xj(y) should be chosen to eliminate the singularity in hj(y). 
An obvious choice is 


, , - h i<y) 

x i (y) = 177T 
h 0 (y) 


(B54) 


The transformation of the independent variable by means of formula (B47) removes 
the singularity from the domain of the problem, because y takes on some small positive 
value when x is zero. Hence, the first-order term hj(y) is nonsingular throughout 
the domain of the problem 0 = x = x^. 
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MATHEMATICAL DEVELOPMENT FOR THE RADIATION-DEPLETED SHOCK LAYER 


The system of equations governing the flow in the stagnation region of a radiating 
shock layer is derived in the text by use of the substitute kernel approximation. This 
system is 


f(T) h'(7)l 


- | eB M (r) - | kpf ( t) h’(r) = 0 


(Cl) 


2f(r?) f"(?7) 



+ a^h(rj) = 0 


(C2) 


f(0) = 0 



= 1 



a 

pxd ~ X) 



(C3) 

(C4) 

(C5) 


h t a) = 1 


(C6) 


This set of equations is subject to the additional condition 


f(r) h’(r) + e 


8 


S A B(t) 
J 0 


where 


r = r w k p 


,C Ta B(t) 
J 0 


3 9 - ? k P T 

dt - I B(r) + | re 2 


-f k P T 

e z dt 


(C7) 


(C8) 


2 

When the parameter e is very much larger than unity and kp, the asymptotic 
solution to the energy equation (Cl) is simply B(r) = Cj + C 2 L Substitution of this solu- 
tion into the asymptotic form of the integral condition (C7) gives C^ = C 2 = 0. This 
solution obviously does not satisfy the boundary condition = 1, which indicates that 

the asymptotic solution is not valid in the vicinity of the boundary (shock) at r = r^. 

This result is not surprising when it is recalled that the existence of a thermal boundary 
layer has been established on physical grounds in the text. 
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j 

f 

ji 

In order to determine the form of the ,T boundary layer" equation valid near the 
{I shock, the stretched coordinate 


/■ 

I 

$ 

j and the functions 

>> 

j 

| 

J; 



B b (4) = B(r) 
%(4) = f(r) 


(C9) 


(CIO) 

(Cll) 


are introduced. The subscript b indicates that these functions are valid only in the 
boundary layer. 

When rewritten in terms of the boundary-layer variables, equation (Cl) becomes 


-e 


3n-l 




+ e 


n-1 9 ,2 


kpf b (|) h b (£) 




If n is set equal to 1 and e is allowed to grow without limit, the most highly 
differentiated term will be retained without losing the significant term of the unstretched 
problem. The resulting differential equation is 


%(?) S> (S >1" +f B b<«> = £ " 2 1 k P f b«> h b<« 


(C12) 


The momentum equation, when written in terms of the stretched coordinate £ = - rjje, 


becomes 


2f b (« f b (?) - fed = -e“ 2 a 2 h b (?) 


(C13) 


It has been shown that the boundary layer is characterized by the parameter e - ^ 
and it would seem proper to expand both the boundary layer and asymptotic solutions as 
power series in e“*. However, f a (r]) (where the subscript a indicates the asymptotic 
solution valid far from the shock) is not analytic in e“* near e - * = 0, but is analytic in 
e -l/2. Consequently, the solutions will be expanded as power series in e - */ 2 ; that is, 


ou 

B a (r) = £ e- n/ 2 B ajI 1 (T) 


n=0 


(C 14) 
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00 


faM = X 6 ^ ^a,n^) 

n=0 

OO 

(C15) 

B b ({) = £ £- n/2 B bi „(4) 
n=0 

OO 

(C16) 

£ b<« = £ «' n/2£ b,n<« 
n=0 

(C17) 


In addition, it will be assumed that the enthalpy h is an analytic function of B 
throughout the interval 0 = B = 1, and from physical considerations it will be assumed 
that h = 0 when B = 0. Then 


h(B) = h(B 0 ) 




2« 

+ -Bjh 



+ . . • 


(C 18) 


where the dot ( ' ) indicates differentiation with respect to the variable Bg. Substitution 
of expansions (C14), (Cl 5), and (C18) into systems (Cl) to (C7) gives 


| Ba,oM + ‘‘ 1/2 | Ba,iM t B^ 2 (t) - |f a>0 M h(B a>0 ) B^ 0 (r) 

+ 1 4 f a,0< T ) h(B a ,o) Ba, 0 ( T | +•••=» 

2 t a,o(>)) 4 a,o(l) - £ a, 0 <l) + a 2 h^B a>0 )> + e” 1/2 { 2f a , 0 ( T )) £ a,l(>l) 

^ L J ^ -v 

2f a,oW f a,l^ + 2f a,0^) f a,lM + &2h ( B a,o) B a,l^} + • • ■ = 0 


f a,0(°) + e“ 1 / 2 f aj i(°) + e -1 f a ,2(0) + • • ■ = 0 


(C19) 


(C20) 

(C21) 
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where 


-fltpt 


r 0 = V^’ 0 B a,0<‘> e 

T 3k t 

r l = r w k pj 0 ’ B a,lW e ^ dt + r w k P T A,l®a,0^ T A,o) 

t _ 3 k t 

r 2 = r w k pJ 0 A ’° B a,2(t) e 2 ? dt + r w k p T Aj2 B aj2 ^ A} ^ + ^kpT^B^r^ 


+ r w k p e 


“2 k P T A,0f er A 


'I 0 a [Bb,o«) - B a,o( T a,o)] d « 


(C23) 

(C24) 

(C25) 

(C26) 

3 ) 

(C27) 


The shock-layer optical thickness is determined from the condition 



APPENDIX C - Continued 


The corresponding boundary-layer equations are obtained by substituting the 
expansions for B b and f b into equations (C12) and (C13) with the result 

f b,0® 1 '( B b,o) B’ bj0 ({) +| B b>0 («)] + a' 1/2 [t bj0 («) b(B b>0 ) 

+ f bj0 <!) h(B b;0 ) B bj0 (l) B b>1 ({) + f b;1 (|) h^o) B^ oW 
* | B b,l<«>]" + 1 '( B b,o) B b,2<» + ° (0*9) 

f 'b,0<« - j * =- 1/2 {2fb,o(?) £,i«> - \ 0 K) ^,1«) 

+ 2f b;0 (O ( b .i(Oj + e-'^ofO f^ 2 (C) + . . j + . . . = 0 (C30) 

f b (0) + e" 1/2 f b}1 (0) + e -1 f bj2 (°) + . • . = o (C31) 



(C35) 
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Systems (C19) to (C24) and (C29) to (C35) lead to a set of recursive systems for 
the solution of B a> i, f a ^, and f^. 


Zero-Order Solutions 

The differential system which describes the asymptotic solutions to zero order in 
the small parameter e“l/2 is 





B a,0^ T ) " 0 

(C36) 

2fa,o0?) 

f a,oM 

- [ f a,0fa)] 2 + a2h ( B a,o) = 0 

(C37) 



f a ,o(0) - 0 

(C38) 


1 

a,O^A,o) = 1 

(C39) 


f a,o( 

A, °^ \j2xd - X) 

(C40) 

B a,0(t) 

-fM 

e 

t-r| „ “| k P T 

dt - 2B a?0 (T) - | r Q e =0 

(C41) 


The solution is 


B a ,C)( T ) = 0 

(C42) 

/ n \ 2 

(C43) 

2 ^ - x) - i * ^<1 -a 

(C44) 

*A,0 

T A,0 _! 

(C45) 


The zero-order boundary-layer system is 

f b,0^) k( B b,o) B b,0^) + f B b,0^) = C l°^ + C 2°^ ( C46 ) 
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2f b,o(« f b,o<0 - [v»] = 0 

(C47) 

f b,0<°> = 1 

(C48) 

lim f b 0 U) = 1 
£ _*.00 ? 

(C49) 

B b>0 (0) - 1 

(C50) 

lim B b q(X) = 0 

I _*oo ’ 

(C51) 

lim B b 0 (£) = 0 

£—°o ’ 

(C52) 

The solution to this system is 


l-ff* «<B)f 
3 jB b,o B 

(C53) 

%,0® - 1 

(C54) 


First-Order Solutions 

The differential system which describes the asymptotic solutions to first order in 
e“ 1/2 i S 

B^j(r) = 0 (C55) 

27?2f a,l M " 4??f a,l^ + 4f a,l( r ?) + 8 xU - x) li(0) B aj i(7?) = 0 (C56) 



(C57) 

(C58) 

(C59) 
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3 

| kpj - ^’ 0 B a>1 (t) e'5 |t ' T|kp dt - 2B M (T) + | rje 2 kpT + § kp^,iB a ,o( T A,o) = 0 


(C60) 


The system for determining B a j(t) is identical to the system for obtaining 
Ba,o( T )’ Therefore, 


B a,l< T > = 0 


(C61) 


and equation (C56) becomes 


2??2f a,l^ " 4??f a,l^ + 4f a,l^ = 0 

subject to the boundary conditions (C57) and (C58). Inspection of the preceding equation 
indicates that condition (C57) is satisfied automatically, so that another independent con- 
dition or equation must be specified in order to obtain a nonarbitrary solution for f aj id). 
This condition can be obtained from the differential system for terms of second order. 

The system which determines the first-order boundary-layer solutions is 

*( B b,o) + h(B bj0 ) B ; j0 d) B M d) + | B bjl d) = C W + cfh (C 62) 


f b>1 (^) = 0 (C63) 

f b>1 (0)=0 (C64) 

lim f b j(^) = 0 (C65) 

£-^°0 9 

B bl (0)=0 (C66) 

lim B b jd) = 0 (C67) 

I— 00 ’ 

lim B: At) = 0 (C68) 

|_*O0 1 


The solution to this trivial system is, of course, 


B b ,id) = 0 


(C69) 
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M«> - 0 


(C70) 


Second-Order Solutions 

The diff erential system which determines the second-order asymptotic solution is 


Ba,2( T ) = 0 


(C71) 



- 

r 1 

2 

47 ? f a,2 (r?) + 4f a,2^) = "^,0 

21a, l(7) f a,l<7> - 

^ — 1 

3 

t-H 

- rt 

+ a 2 h(0) B a>2 (r7) 


f a f 2(0) = 0 


l, 2 ^A,o) - ^A,l f a, 1(^,0) 


? A,2 I' A, 1 


^A,Q j 1^,0 1 


f a,2^A 5 o) " ^A,l f a, l(^A,'o) " 2 ^ 


A, 2 


(C72) 

(C73) 

(C74) 

(C75) 


| kp^ A ’° B a 2 (t) e 2 kp 4 dt - 2B a>2 (r) + | r 2 e 

f k p( T A,(T 


A,0 

— k r 
2 p 3 


+ f k P T A,l B a,l( T A,o) 


+ 4 k P r A 9! B a.o( T A.n) + 9 k p e 


2 ft P A,2 a,0^ A,0j 
The solution to this system is 

B a ,2( T ) = 


I 0 A B b,0(^) - B a,o( T A ,o) 


d£ = 0 (C76) 


,(i + r w ) 


l+-kpT 


2 1 +|(l - r w) k P T A,0 


(C77) 


The second condition for the quantity f & 1 (t?) can be obtained by evaluating equa- 
tion (C72) at 77 = 0, which gives 


E^^O) = aJli(O) B a?2 (0) 


(C78) 
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With this condition, the solution for f aj i (??) can be completely specified with the result 




In addition, 


and 


%1 = - a 4,0^(°) B a,2<°) 
%1 


(C79) 


(C80) 


Kp(0) 

As before, it is necessary to look to the next higher order system in order to 
obtain a second independent condition or equation for f a 2 (^). This condition is 


(C81) 


f M<°> » 0 


(C82) 


In order to solve equation (C72) for f a 2(v)> it is necessary to express the optical 
path length r as a function of rj. This relationship can be made with aid of the defi- 
nitions of r and T], The result is 


t = t a n —2— + o(e-l) 
A ’° «A,0 V ’ 

Now, equation (C72) can be written in the form 

2 2 2f a,2^) _ 4r ? f a,2^) + 4f a,2^) = -a 2 r7 AjQ h(0)- 


(C83) 


3 (l + r w) k p T A,0 


1 + 


!C - r ») 


k P r A,0 


The solution to this equation is 

,2 


f a,2W = 


V 


iV 


f A,0 j 


V A,2 

. 2 _. A loge „ i;0 


ifl 


A l0 Se ^AjO 


V 


V + A — ^ — log. 


V 


'a,0 ^,0 


'A,0 


(C84) 
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where 


Also, 


and 




2 

) k P T A,0 


Ni 

^A,2 - -1A,( fe 


2r 

) 1 + f k P T A,0 ^ “ lo ^e ^A,oj 


\2 

r A,2 = Tw + ^A,0 

’ Kp (0) 


^(0) 


^ t *,0 +€1/2t A,1 +T A,2| 


fA,( 


N b >, 2 (°>Hva,«) 

k l (o> v 7 


K p X ( B b,0^j) _ 1 


/Cp(O) 




(C85) 


(C86) 


Radiative Flux and Standoff Distance 


In the section "The Optically Thick Shock Layer," it was shown that with the sub- 
stitute kernel approximation, the radiant flux passing into the cold wall can be expressed 
in terms of the black-body emissive power of the gas adjacent to the wall. Substituting 
the solution for B a (0) into the expression (171) gives 


R 

q w 


1 - r. 


w 


1 + 


f (‘ - r w) 


k P T A,0 


(C87) 


The ratio of the shock standoff distance to the shock standoff distance for radiation- 
less flow is given by the condition 


A = 



Kp(B) dt 


(C88) 
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Expanding this expression in powers of gives 


A — It' 


r 

• 



^A 0 

h(0) + h(0p-- — L 

V- 

/cp(0) 


r eT A,0^ 1//2r A,l +T A,2 
Jn 


B a ,2<0)(l + |kpT a _ 0 j 


h ( B b,0<a) *p ' ( B b,0<«) d « 


(C89) 
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APPENDIX D 
SYMBOLS 

a constant defined by equation (55) 

a n nth-order coefficient in perturbation expansion of constant a 

a* constant defined by equation (145) 

B black-body emissive power, nondimensional except in section "Stagnation 

Model for a Radiating Shock Layer" where it is 2 T^,' erg/cm^-ster-sec 

B^ Planck function, nondimensional except in section "Stagnation Model for a 

Radiating Shock Layer" where it is defined by equation (11), 
erg/cm3 -ster-sec 

B nondimensional black-body emissive power in section "Stagnation Model for a 

Radiating Shock Layer" 

B^ nondimensional Planck function in section "Stagnation Model for a Radiating 

Shock Layer" 

B 0 O=B[h 0 ()] 

Bx,0( ) - B X [h 0 ( )] 

B 2 value of nondimensional black-body emissive power in interior of an optically 

thick shock layer 

B w value of nondimensional black-body emissive power in gas adjacent to wall in 

an optically thick shock layer 

> 

B,B . constants defined in section "The Optically Thick Shock Layer" 

B a nondimensional black-body emissive power in interior of a radiation-depleted 

shock layer 

B a?n nth-order coefficient in perturbation expansion of B a 
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nondimensional black-body emissive power in shock boundary layer in a 
radiation-depleted shock layer 

nth-order coefficient in perturbation expansion of 

constant in formula (240) 

constant defined by equation (193) 

constant defined by equation (149) 


''n 


?> 


constants of integration 


j 


velocity of light, cm/sec 

exponential integral function of order n 

functions of h defined by equation (55) 
nondimensional stream function defined by equation (36a) 
nth-order coefficient in perturbation expansion of f 
nth-order coefficient in P-L-K expansion of f 

nondimensional stream function in interior of a radiation-depleted shock 

nth-order coefficient in perturbation expansion of f a 

nondimensional stream function in shock boundary layer in a radiation- 
depleted shock layer 

nth-order coefficient in perturbation expansion of f^ 

stream function; used in section "Stagnation Model for a Radiating Shock 
Layer,” g/cm^-sec 



APPENDIX D — Continued 


nondimensional stream function in viscous boundary layer; defined by 
equation (86) 

nth-order coefficient in perturbation expansion of g 

enthalpy, nondimensional except in section "Stagnation Model for a Radiating 
Shock Layer" where it is the static specific enthalpy, erg/g; also Planck's 
constant, erg/sec 

nth-order coefficient in perturbation expansion of h, the nondimensional 
enthalpy 

nth-order coefficient in P-L-K expansion of h, nondimensional enthalpy 

nth-order coefficient in expansion of h, the static specific enthalpy, erg/g 

average value of nondimensional enthalpy h (eq. (135)); also nondimensional 
enthalpy in section "Stagnation Model for a Radiating Shock Layer" 

nth-order coefficient in perturbation expansion of h 

value of h for which y ^exponent in correlation formula /Cp(h) = Ch^j 
changes 

total enthalpy in section "Stagnation Model for a Radiating Shock Layer," 

erg/g 

value of nondimensional enthalpy in interior of an optically thick shock layer 

enthalpy in gas adjacent to wall, nondimensional except in section "Stagnation 
Model for a Radiating Shock Layer" where it is the static specific enthalpy 
in gas at wall conditions, erg/g 

value of nondimensional enthalpy at outer edge of viscous boundary layer 

divergence of radiant flux vector, nondimensional except in section "Stagnation 
Model for a Radiating Shock Layer" where it has the dimensions, 
erg/cm^-sec 
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nth-order coefficient in perturbation expansion of nondimensional divergence 
of radiant flux vector, I 

nondimensional enthalpy in viscous boundary layer; defined by equation (85) 
nth-order coefficient in perturbation expansion of i 

nondimensional divergence of radiant flux vector in viscous boundary layer 

nth-order coefficient in the perturbation expansion of J 

specific intensity of radiation, erg/cm3-ster-sec 

mass emission coefficient, erg/g-cm-ster-sec 

Boltzmann's constant, erg/°K 

effective coefficient of heat conduction including energy transport by 

molecular collisions and by diffusion of reacting species, erg/cm-sec-°K 

Bouguer number P s K p g A A 

quantity defined in equation (107) 

direction cosine between direction of a beam of intensity and 
ith direction 

Boltzmann number 

Peclet number 

Prandtl number 

Reynolds number 

pressure, dyne/cm 2 

standard pressure of air, 1.013 x 10 6 dyne/cm 2 
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(n) 


u. 

qR 

q? 

R 

V 

1w 

R 

% 

R s 


r 0,X 


1,X 


L w 


nth-order coefficient in expansion p, dyne/cm2 

ith component of combined radiant and conductive heat fluxes, erg/cm2-sec 

nth-order coefficient in expansion of q^, erg/cm2-sec 

component of conductive heat flux vector in ^-direction, erg/cm2-sec 

component of conductive heat flux vector in ith direction, erg/cm2-sec 

component of radiant heat flux vector in 77-direction, erg/cm2-sec 

component of radiant heat flux vector in ith direction, erg/cm2-sec 

component of monochromatic radiant heat flux vector in ith direction, 
erg/cm2-sec 

nondimensional rate of radiant heat transfer to wall 
gas constant for air, 2.882 x 10®, cm2/sec2-°K 
body nose radius, cm 

shock radius in vicinity of stagnation streamline, cm 

position coordinate, cm 

defined by equation (108) 

defined by equation (113) 

reflectivity of wall 

radiation source function, erg/em3-ster-sec 
position coordinate, cm 
nondimensional variable of integration 
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x 0 

y 

z 

a 

P 


APPENDIX D - Continued 

temperature, °K 

standard temperature, 273.16° K 

temperature of wall, °K 

dummy variable of integration 

nondimensional variable of integration 

component of gas velocity in r-direction, cm/sec 

nth-order coefficient in expansion of u, cm/sec 

component of gas velocity in ith direction, cm/sec 

volume, cm3 

free-stream velocity, cm/sec 

component of gas velocity in z-direction, cm/sec 

nth-order coefficient in expansion of w, cm/sec 

coordinate in transformed plane; also normalized Dorodnitsyn coordinate, 

Vvo 

nth-order coefficient of P-L-K expansion of x 
value of transformed coordinate for which 77 = 0 
coordinate in transformed plane 
position coordinate, cm 

magnitude of step in step function model of absorption coefficient 
local angle of inclination of bow shock from stream direction 
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y 

y 


W2 

A 

A A 

A 


6 


mass extinction coefficient, cm2/g 

Euler’s constant, y = 0.577216. . . 

exponent in correlation formula, K-p = 

exponents in correlation formula for /Cp (eq. (148)) 

shock standoff distance, cm 

shock standoff distance for nonradiating shock layer, cm 

ratio of shock standoff distance for radiating and nonradiating shock layer, 
A/A A 

displacement distance for viscous boundary layer, cm; also exponent in cor- 
relation formula B = h® 


5 n nth-order coefficient in perturbation expansion of displacement distance, 

6 , cm 

4/3 

e radiation cooling parameter, 4aT g kp/p oo W 00 

£ transformed nondimensional Dorodnitsyn coordinate in radiation-depleted 

shock layer 

£ transformed optical path length in optically thick shock layer 

77 Dorodnitsyn coordinate, nondimensional except in section ’’Stagnation Model 

for a Radiating Shock Layer’’ where it is defined by equation (35), g/cm^ 

77 nondimensional Dorodnitsyn coordinate in section "Stagnation Model for a 

Radiating Shock Layer" 

77 ’ variable of integration 

77 ^ location of shock in terms of Dorodnitsyn coordinate, nondimensional except 

in section "Stagnation Model for a Radiating Shock Layer" where it has 
dimensions, g/cm 2 
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nondimensional location of shock in terms of Dorodnitsyn coordinate in sec- 
tion "Stagnation Model for a Radiating Shock Layer" 

nth-order coefficient in perturbation expansion of 77^ 

nth-order coefficient in P-L-K expansion of 17 

angle defined in sketch (a) 

mass absorption coefficient, nondimensional except in section "Stagnation 
Model for a Radiating Shock Layer" where it has dimensions, cm^/g 

nondimensional mass absorption coefficient in section "Stagnation Model for 
a Radiating Shock Layer" 

Planck mean mass absorption coefficient, nondimensional except in section 
"Stagnation Model for a Radiating Shock Layer" where it has dimensions, 
cm2/g 

nondimensional Planck mean mass absorption coefficient in section "Stagna- 
tion Model for a Radiating Shock Layer" 

evaluated at h = 1 


nth-order coefficient in perturbation expansion of Kp 
Rosseland mean mass absorption coefficient, cm^/g 

■" 1/2 

wavelength, cm; also boundary-layer parameter, ^Np e ^ 
coefficient of viscosity, dyne-sec/cm2 
second coefficient of viscosity, dyne-sec/cm2 

transformed nondimensional Dorodnitsyn coordinate in viscous boundary layer 


thickness of viscous boundary layer in terms of £ 
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P 


p o 

p 2 

p( n ) 


a 



T 


T 


X 



T 


A 


T A,n 

T A,A 
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density, g/cm^ 

standard density, 1.288 x 10"^ g/cm^ 

density in interior of an optically thick shock layer, g/cm3 
nth-order coefficient in expansion of p, g/cm^ 

Stefan-Boltzmann constant, 5.669 x 10“^ erg/cm.2-sec-°K^ 
area, cm^ 

mass scattering coefficient, cm2/g 

transformed monochromatic optical path length in viscous boundary layer 
thickness of viscous boundary layer in terms of 

optical path length in a gray gas, normalized except in section "Stagnation 
Model for a Radiating Shock Layer" 

normalized optical path length in a gray gas in section "Stagnation Model for 
a Radiating Shock Layer" 

monochromatic optical path length, normalized except in section "Stagnation 
Model for a Radiating Shock Layer" 

normalized monochromatic optical path length in section "Stagnation Model 
for a Radiating Shock Layer" 

monochromatic optical path length in s-direction 

thickness of viscous boundary layer in terms of 

shock location in terms of r 

nth-order coefficient in perturbation expansion of 
shock location in terms of r. 
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X 
co 

o>i 
W 2 

Subscripts: 

b shock boundary layer in radiation-depleted solution 

s indicates value of dimensional quantity at normal-shock equilibrium 

conditions 

00 indicates value of dimensional quantity in the free stream 

0,1 order of perturbation 

Primed symbols denote differentiation with respect to the argument. An arrow 
over a symbol denotes a vector. Bars over symbols denote nondimensional quantities in 
section "Stagnation Model for a Radiating Shock Layer." 


component of viscous stress tensor, dyne/cm^ 
functions of 77 defined by equations (109) and (114) 
density ratio across normal shock, Pao/pQ 
solid angle 

constant defined by equation (188) 
constant defined by equation (197) 


% °> 




I . 


X 
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Figure 2.- Temperature correlation. T = ^p/ P q) °' 09 f 2 (h/RT q) , °K. 







(a) p/p 0 = Iff 3 . 

Figure 3.- Planck mean mass absorption coefficient correlation. 


122 


n- 





t 


0 Kivel and Pailey (ref. 43) 

□ Nardone et al. (ref. 50) 

O Sewell (ref. 46) 

A Gilmore (ref. 44) 

IA Armstrong et al. (ref. 45) 


i 

l 




T, °K 



(b) p/p o “ 10 -2 . 
Figure 3.- Continued. 
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(c) p/p 0 = 10’ 1 . 
Figure 3.- Continued. 


124 





<d> p/p 0 = 10°. 
Figure 3.- Continued. 
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(e) p/p 0 = 10 1 . 
Figure 3.- Concluded. 
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Figure 4.- Comparison of the Poincare-Lighthi I l-Kuo (P-L-K), conventional perturbation, and exact solutions, 
e = 0.1, kp = 0, p = Constant, KpB = h<>. 
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Figure 9.- Effect of the enthalpy variation of the absorption coefficient on the rate of radiant heat transfer 

to the stagnation point, e = 0.1; r w = 0. 
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Spectral distribution Cref. 50) 



X, micrometers 

Figure 13.- Step function model of the mass absorption coefficient of high temperature air. 
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Figure 16.- Comparison of the optically thin and small perturbation solutions with numerical results of reference 9. 

W„ = 9.75 km/sec; p s = 10 atm-, R s = 1.525 m. 
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Trajectory 1 Iron sphere-vertical entry at 70 km/sec, R £ 
Trajectory 2 Iron sphere-vertical entry at 70 km/sec, R £ 
Trajectory 3 Mars probe-shallow entry at 15 km/sec, R g 
Trajectory 4 Mars probe-shallow entry at 15 km/sec, R s 


Figure 31.- Entry trajectories in the £,kp space. 
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